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The  characteristics  of  the  coherent  structures  of  a  plane  jet 
were  evaluated  experimentally  and  with  a  kinematic  numerical 
representation . 
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frequency  were  measured  for  10<x/D£60.  These  properties,  including 
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do  ,:sr„  1473  EDITION  OF  1  NOV  SS  It  OBSOLETE 

S/N  0103-014* t*Sl  I  ..  — — — g-  .  R'wnm.  ■*  wr 


SECURITY  CLASSIFICATION  OF  THIS  R« 


ifcJi 


self-preserving  flow  variables. 

Spatially  coherent  patterns  within  the  turbulent  velocity 
field  were  evaluated  utilizing  two  point  velocity  correlation 
measurements.  From  these,  it  was  found  that  the  instantaneous 
velocity  fluctuation  distributions  are  antisymmetric  with  respect 
to  the  centerline  and  coherent  over  the  entire  width  of  the  flow 
field.  In  addition,  the  patterns  are  periodic  in  the  streamwise 
direction. 

The  two-dimensionality  of  the  large  structures  was  evaluated 
through  the  comparison  of  longitudinal  and  vertical  integral  lengtl 
scales  derived  from  correlations  in  both  the  fully  turbulent  zone 
and  the  intermittent  region.  In  the  intermittent  region,  the 
integral  scales  were  approximately  equal  while  in  the  fully 
turbulent  zone,  the  longitudinal  scale  was  approximately  double 
the  vertical  one. 

Utilizing  a  combination  of  flow  visualization,  structural 
passage  frequency  measurement  and  similarity  considerations,  the 
convective  velocity  of  the  vortex-like  structures  was  estimated  to 
be  one  half  the  mean  centerline  velocity. 

From  the  experimental  results,  a  hypothetical  vortex  street 
representation  of  the  plane  jet  was  formulated  and  evaluated 
numerically.  This  kinematic  representation  was  found  to  be 
capable  of  generating  realistic  distributions  of  mean  longitudinal 
and  lateral  velocity  as  well  as  the  proper  spreading  and  velocity 
decay  rates.  Although  the  correct  magnitude  of  the  Reynolds 
stress  was  calculated,  its  sign  was  in  error.  The  small  scales 
of  turbulence  were  not  included  in  the  model  and  consequently, 
it  was  not  possible  to  predict  reasonable  distributions  of  the 
fluctuation  intensities,  the  energy  spectra  or  interface 
properties . 
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Qler,  James  Walter,  Ph.D.,  Purdue  University,  May  1980. 
Coherent  Structures  in  the  Similarity  Region  of  a  Two- 
dimensional  Turbulent  Jet  :  A  Vortex  Street.  Major 
Professor:  V.  W.  Goldschmidt. 

^  The  characteristics  of  the  coherent  structures  of  a 
plane  jet  were  evaluated  experimentally  and  with  a  kinematic 
numerical  representation. 

Distributions  of  intermitWency  /and  interface  crossing 


frequency  were  measured  for  10*x/D^60.  These  properties, 
including  the  maximum  interface  crossing  frequency,  were 
found  to  scale  as  self-preserving  flow  variables. 

Spatially  coherent  patterns  within  the  turbulent 
velocity  field  were  evaluated  utilizing  two  point  velocity 
correlation  measurements.  From  these,  it  was  found  that 
the  instantaneous  velocity  fluctuation  distributions  are 
antisymmetric  with  respect  to  the  centerline  and  coherent 
over  the  entire  width  of  the  flow  field.  In  addition,  the 
patterns  are  periodic  in  the  streamwise  direction. 

The  two-dimensionality  of  the  large  structures  was 
evaluated  through  the  comparison  of  longitudinal  and 
vertical  integral  length  scales  derived  from  correlations 
in  both  the  fully  turbulent  zone  and  the  intermittent 
region.  In  the  intermittent  region,  the  integral  scales 
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were  approximately  equal  while  in  the  fully  turbulent  zone, 
the  longitudinal  scale  was  approximately  double  the  vertical 
one. 

Utilizing  a  combination  of  flow  visualization,  struc¬ 
tural  passage  frequency  measurement  and  similarity  consid¬ 
erations,  the  convective  velocity  of  the  vortex-like 
structures  was  estimated  to  be  one  half  the  mean  centerline 
velocity. 

From  the  experimental  results,  a  hypothetical  vortex 
street  representation  of  the  plane  jet  was  formulated  and 
evaluated  numerically.  This  kinematic  representation  was 
found  to  be  capable  of  generating  realistic  distributions 
of  mean  longitudinal  and  lateral  velocity  as  well  as  the 
proper  spreading  and  velocity  decay  rates.  Although  the 
correct  magnitude  of  the  Reynolds  stress  was  calculated, 
its  sign  was  in  error.  The  small  scales  of  turbulence  were 
not  included  in  the  model  and  consequently,  it  was  not 
possible  to  predict  reasonable  distributions  of  the  fluc¬ 
tuation  intensities,  the  energy  spectra  or  interface 
properties. 


CHAPTER  I 


INTRODUCTION 


In  recent  years,  there  has  been  a  gradual  shift  in 
the  approach  to  turbulent  fluid  mechanics  research.  It  is 
becoming  increasingly  clear  that  much  of  what  has  been 
traditionally  considered  random  phenomena  has  much  more 
deterministic  aspects  than  previously  thought.  In 
particular,  a  growing  body  of  experimental  evidence  is 
accumulating  which  demonstrates  that  much  of  the  physics  of 
turbulent  shear  flow  may  be  controlled  by  a  finite  set  of 
coherent  large  scale  eddying  motions.  The  relati.e 
importance  of  these  large  scales  has  been  emphasized  in 
several  recent  review  articles  (Kovasznay  (1977) ,  Roshko 
(1976),  Davies  and  Yule  (1975),  Laufer  (1975)  and  Mollo- 
Christensen  (1971) )  .  While  the  acceptance  of  the  nonrandom 
aspects  of  turbulence  is  common,  it  is  by  no  means  universal. 
(See  for  instance  Pui  and  Gartshore  (1979) ,  Chandrsuda,  et 
al.  (1978)  and  Bradshaw  (1976)).  However,  the  presence  of 
the  large  scale  structures  and  to  some  extent,  their 
importance,  have  become  experimentally  established  facts. 

The  goal  of  the  research  described  herein  has  been  to 
examine  the  fully  developed  two-dimensional  plane  turbulent 
jet  for  evidence  of  the  effects  of  the  large  scale  eddies 
characteristic  of  that  flow  field.  One  of  the  original 
incentives  was  to  investigate  the  transport  characteristics 
within  the  jet  conditioned  on  the  instantaneous  position  of 
the  turbulent-nonturbulent  interface.  During  the  course  of 
these  measurements,  the  presence  of  an  apparent  antisymmetric 


and  periodic  vortical  pattern  was  noted,  suggesting  a  two- 
dimensional  and  ordered  pattern  similar  to  that  found  in 
two-dimensional  mixing  layers  (Brown  &  Roshko  (1971))  and 
wakes  (Papailiou  and  Lykoudis  (1974)).  Consequently,  the 
majority  of  the  experiments  and  modeling  effort  was  devoted 
to  evaluating  the  degree  to  which  these  features  are  present 
in  the  plane  turbulent  jet. 
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CHAPTER  II 
GENERAL  BACKGROUND 


The  purpose  of  this  chapter  is  to  first  introduce  the 
reader  to  the  classical  characteristics  and  descriptive 
terms  for  turbulent  plane  jets.  Following  that,  a  survey  of 
the  literature  pertaining  to  large  scale  structures  in  other 
free  turbulent  shear  flows  will  be  given.  Finally,  the 
concept  of  a  von  Kcirmcin  'vortex  street'  type  pattern  will  be 
developed  for  a  two-dimensional  plane  turbulent  jet. 

A.  The  Turbulent  Plane  Jet 
1.  Initial  Region 

Turbulent  plane  jets  result  from  the  flow  of  a  two- 
dimensional  nozzle  emanating  into  an  infinite  environment. 

The  most  general  case  is  that  resulting  from  an  initially 
uniform  velocity  profile  (except  for  two  boundary  layers) 
flowing  into  an  otherwise  quiescent  environment.  This  is  the 
case  now  treated. 

A  plan  view  of  the  flow  field  is  given  in  Figure  2-1. 
Three  distinct  regions  of  flow  development  may  be  defined, 
the  first  of  which  is  the  initial  region  beginning  at  the 
jet  exit.  The  initial  region  is  characterized  by  a  potential 
cone  of  uniform  velocity,  UQ,  bounded  on  either  side  by 
plane  turbulent  shear  layers. 

For  properly  designed  nozzles,  the  shear  layers 
originate  from  initially  laminar  boundary  layers,  then 
transition  to  turbulence  and  develop  linearly  with  respect 
to  downstream  distance.  The  growing  mixing  layers  eventually 
meet  at  the  centerline  which  marks  the  end  of  the  potential 
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cone  and  initial  region.  This  usually  occurs  four  to  six 
slotwidths,  D,  from  the  jet  exit.  Prior  to  this  point,  the 
shear  layers  interact  through  the  mutually  induced  pressure 
field  only. 

2 .  Transition  Region 

From  the  point  of  contact  at  the  end  of  the  potential 
cone,  the  shear  layers  go  through  a  period  of  interaction  or 
transition  until  a  state  of  dynamic  equilibrium  is  reached 
across  the  entire  width  of  the  flow  field.  This  is 
generally  accomplished  within  the  first  30  slot  widths  of 
flow  travel. 

3.  Similarity  Region 

Similarity  exists  for  a  developing  flow  field  when 

one  velocity  scale  and  one  length  scale  are  sufficient  to 

make  distributions  of  reduced  velocities  in  terms  of 

reduced  coordinates  identical.  The  two-dimensional  plane 

jet  is  characterized  by  this  property  as  demonstrated  for 

distributions  of  the  mean  longitudinal  velocity,  U,  in 

Figure  2-2.  The  velocity  and  length  scales  utilized  are  the 

mean  centerline  velocity,  and  the  velocity  halfwidth,  b 

(the  distance  from  the  jet  centerline  to  the  lateral 

location  at  which  U  =  .5U  ).  From  derivations  based  on 

m 

similarity  arguments  (See,  for  instance,  Townsend  (1956)  or 

Hinze  (1975)),  it  may  be  shown  that  the  velocity  and  length 

scales  must  themselves  be  functions  of  the  downstream 
. .  -1/2 
distance,  x.  In  particular,  one  expects  Um  «  (x/D)  and 

b  *  (x/D)  which  is  verified  by  the  data  illustrated  in 

Figures  2-3  and  2-4.  The  variation  of  these  scales  is 

independent  of  the  jet  Reynolds  number,  =  UqD/v.* 


*  In  a  strict  sense,  perfect  similarity  requires  an  inde¬ 
pendence  of  the  initial  conditions  which  is  only  approxi¬ 
mately  true  in  real  flows.  This  is  evidenced  by  a  lack 
of  coincidence  of  the  kinematic  and  dynamic  virtual 
origins. 


Figure  2-3  0 /U  vs.  x/D 
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One  of  the  most  distinctive  characteristics  of 
turbulent  shear  flows  with  at  least  one  free  boundary  is  the 
existence  of  a  sharply  defined  interface  separating  the 
vortical,  or  turbulent  flow,  from  the  nonvortical,  or  non- 
turbulent,  ambient  fluid.  Driven  by  the  large  scales  of 
the  turbulent  motion  (Corrsin  &  Kistler  (1954)),  the  inter¬ 
face  is  characterized  by  large  bulges  whose  maximum  and 
minimum  extensions  define  the  width  of  the  intermittent 
region.  For  an  observer  at  a  fixed  point  in  the  intermit¬ 
tent  region,  the  flow  may  be  described  as  intermittently 
turbulent  as  the  deformations  in  the  turbulent-nonturbulent 
interface  are  convected  past.  Following  the  convention  of 
Townsend  (1949),  an  intermittency  function,  I(r,t),  may  be 
defined  as  taking  on  the  values  of  one  or  zero  when  the  flow 
is  turbulent  or  nonturbulent ,  respectively  at  a  fixed  point, 
r,  in  the  flow  field.  A  temporal  average  of  I(r,t)  is  the 
intermittency,  y(r),  or  fraction  of  time  that  the  flow  is 
turbulent  at  r.  The  interface  crossing  frequency,  f^,,  is 
defined  as  the  rate  at  which  the  bulges  in  the  interface 
cross  a  fixed  point. 

B.  Large  Scale  Structures  in  Turbulent 
Free  Shear  Flows 

Recognition  of  the  importance  of  the  large  scale 
structures  in  turbulent  shear  flows  dates  back  at  least  as 
far  as  Townsend  (1956)  and  Grant  (1958).  Townsend's  concept 
of  a  double  structured  nature  of  turbulence  and  hypothesis 
for  the  "dynamic  equilibrium  of  the  large  structures"  remain 
quite  useful  in  the  interpretation  of  recent  data  concerning 
the  large  structures.  He  postulated  that  the  large  scale 
motion  plays  a  central  role  in  the  dynamics  of  flow  develop¬ 
ment  and  energy  exchange  between  the  mean  flow  and  small 
scales.  A  state  of  equilibrium  exists  between  the  rate  at 
which  energy  is  extracted  from  the  mean  flow  by  the  large 
scale  structures  and  the  rate  at  which  these  structures  lose 


energy  to  the  small  scales  of  turbulence. 

Coherent  turbulent  structures  have  since  been 
identified  in  both  bounded  and  free  shear  flows  although 
the  boundary  conditions  and  instability  mechanisms  are 
different. 

Recent  revelations  concerning  the  coherent  motions 
near  the  wall  of  a  boundary  layer  have  somewhat  confused 
the  concept  of  a  "cascade"  of  energy  from  the  large  to  small 
scales  since  the  relationship  appears  to  be  reversed  for 
that  flow.  However,  in  free  shear  flows,  the  energy 
cascade  and  dynamic  equilibrium  hypothesis  appear  to  be 
valid.  The  following  discussion  will  deal  only  with  large 
scale  structures  as  they  relate  to  free  shear  flows,  such  as 
the  plane  turbulent  jet. 

1.  Turbulent  Mixing  Lavers 

Much  of  the  current  emphasis  on  the  large  eddies  of 
turbulence  is  a  result  of  the  flow  visualization  experiments 
of  Brown  &  Roshko  (1971,  1974)  in  a  turbulent  mixing  layer. 
Their  shadowgraphs  of  the  mixing  layer  between  two  parallel 
streams  of  unequal  velocity  and  density  clearly  illustrate 
the  existence  of  large,  essentially  two-dimensional  roller 
vortices  with  axis  perpendicular  to  the  plane  of  mean  shear. 
Photographs  indicate  that  the  large  structures  scale  on  the 
dimensions  of  the  mean  flow  as  they  are  convected  downstream 
at  a  velocity  approximately  equal  to  the  average  of  the  two 
stream  velocities.  In  fact,  it  is  apparent  that  the  growth 
of  the  mixing  layer  is  controlled  by  the  large  eddies  as 
they  undergo  discontinuous  development  through  a  mechanism 
of  eddy  coalescence.* 


*  The  coalescence  event  has  been  alternately  described  as 
merger,  pairing  and  aglomeration.  The  term  "coalescence" 
will  be  used  throughout  this  work. 


Visualization  experiments  were  repeated  over  a  wide 
range  of  Reynolds  numbers  with  no  apparent  effect  on  the 
growth  rate  or  eddy  coalescence.  The  only  difference  in 
the  flows  was  the  increasingly  high  small  scale  turbulent 
energy  for  higher  Reynolds  numbers.  This  is  explained  by 
Townsend's  dynamic  equilibrium  hypothesis  in  that  as  more 
energy  is  extracted  by  the  large  scales  from  the  mean  flow 
for  the  higher  Reynolds  numbers,  an  equally  greater  amount 
of  energy  is  transferred  from  the  large  eddies  to  the  smaller 
ones.  This,  in  turn,  requires  increasingly  smaller  scales 
for  the  viscous  dissipation  of  energy  as  described  by  the 
Kolmogoroff  hypothesis.  Hence,  the  large  eddies  are  inde¬ 
pendent  of  viscosity  and  dependent  only  on  the  mean  flow 
characteristics .  The  small  scales,  on  the  other  hand,  are 
dependent  on  viscosity  due  to  the  overall  energy  dissipation 
requirements. 

Brown  &  Roshko  emphasized  that  their  mixing  layer  fit 
the  classical  characteristics  of  fully  developed  turbulent 
mixing  layers  in  every  way,  including  mean  and  fluctuating 
velocities  and  stress  distributions.  The  large  eddies  are 
believed  to  be  characteristic  of  the  flow  rather  than  the 
result  of  some  artificial  experimental  forcing. 

Winant  &  Browand  (1974)  and  Browand  &  Weidman  (1976) 
made  similar  observations  for  a  mixing  layer  in  water  at  low 
Reynolds  numbers.  In  the  latter  investigation,  a  combination 
of  flow  visualization  and  conditional  sampling  was  used  to 
examine  the  details  of  the  vortex  coalescence  process.  By 
ensemble  averaging  the  velocity  signals  resulting  from  the 
passage  of  pairs  of  vortices  in  various  transitional  states 
of  coalescence,  the  instantaneous  characteristics  of  the 
process  were  deduced.  It  was  found  that  the  vorticity  of  an 
isolated  noninteracting  vortex  is  concentrated  near  its 
center  and  that  the  same  vorticity  distribution  is  recovered 
following  the  coalescence  of  two  such  isolated  vortices.  In 
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addition,  the  rate  of  Reynolds  stress  production  during  the 
active  coalesence  of  vortices  was  found  to  be  almost  double 
that  for  a  passive  vortex. 

An  interesting  argument  was  developed  by  Browand  6 
Weidman  to  emphasize  the  importance  of  the  large  eddies  to 
Reynolds  stress  production.  They  noted  that  their 
distribution  of  measured  dimensionless  Reynolds  stress  is 
essentially  identical  to  that  of  Liepmann  &  Laufer  (1947) 
in  spite  of  the  large  difference  of  Reynolds  number.  This 
is  not  true  for  the  distributions  of  the  r.m.s.  velocity 
fluctuations.  Based  on  these  observations,  they  suggested 
that  it  must  be  the  large  scale  eddies  which  are  responsible 
for  the  bulk  of  the  Reynolds  stress  production  since  it  is 
only  those  which  are  unaffected  by  the  Reynolds  number. 

In  response  to  skepticism  concerning  the  universality 
of  the  two-dimensional  structure  in  mixing  layers  with 
Reynolds  numbers  (see  for  instance,  Bradshaw ( 1976) ) , 

Dimotakis  &  Brown  (1976)  performed  an  investigation  of  this 
effect.  While  Brown  &  Roshko  utilized  a  mixing  layer  at 

5 

Re  =  10  and  Wmant  &  Browand  used  Re  =  300,  the  mixing 
layer  of  Dimotakis  &  Brown  was  investigated  at  Re  *  3  x  10®. 
The  two-dimensional  structure  and  vortex  coalescence 
mechanism  were  noted  to  still  be  apparent. 

Dimotakis  &  Brown  also  investigated  the  entrainment 
mechanism  by  observing  the  paths  of  dye  streaks  injected 
into  the  free  streams.  They  found  that  entrainment  is 
accomplished  in  two  distinct  stages; 

1)  Dye  is  drawn  into  the  mixing  layer  symmetrically  from 
both  sides  through  an  engulfment  process  induced  by  the  large 
structures. 

2)  Once  drawn  into  the  large  eddies,  little  molecular 
mixing  is  evident  since  the  entwined  dye  streaks  remain 
identifiable.  This  situation  perseveres  until  the  identity 
of  the  eddy  is  lost  through  coalescence  or  a  secondary 
mechanism  they  described  as  vortex  tearing.  At  that  point 
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complete  diffusion  of  the  dye  is  almost  instantaneous. 

The  concept  of  an  essentially  two-dimensional  structure 
in  a  fully  developed  turbulent  flow  has  not  met  with 
universal  acceptance.  Chandrsuda,  et  al.  (1978)  and  Pui  6 
Gartshore  (1979)  reported  experiments  in  mixing  layers 
where  the  large  scale  two-dimensional  eddies  either  were  not 
evident  at  all  or  transitioned  to  a  fully  three-dimensional 
turbulent  field  a  short  distance  from  the  splitter  plate. 
Degradation  of  the  two-dimensional  structure  appeared  to  be 
caused  by  comparatively  low  level  turbulence  in  the  free 
streams  or  in  the  initial  shear  layer.  The  proponents  of 
"organized"  turbulence  were  undaunted,  however,  and 
Wygnanski,  et  al.  (1979)  reported  an  equally  well  documented 
experiment  in  which  the  two-dimensional  structure  survived 
severe  background  turbulence  in  the  free  streams  and  in 
the  initial  shear  layer. 

2.  Wake  Flows 

The  existence  of  a  well  organized  vortex  street  in  the 
near  wake  of  a  cylinder  has  long  been  known,  but  until 
recently,  this  structure  was  thought  to  degenerate  quickly 
with  downstream  distance  (see  Roshko  (1954)).  Experiments 
by  Taneda  (1959),  Bevilaqua  &  Lykoudis  (1971),  Papailiou  6 
Lykoudis  (1974)  and  Townsend  (1979)  have  indicated  that  the 
large  eddies  of  the  turbulent  motion  may  be  characterized  by 
the  von  Ketrm^n  vortex  street  at  much  greater  downstream 
distances  and  Reynolds  number  than  previously  thought 
possible. 

Although  the  vortex  street  in  a  plane  wake  is  a 
natural  intuitive  extension  of  the  roller  vortices  in  a 
mixing  layer,  there  appear  to  be  fundamental  differences. 
First,  in  none  of  the  experiments  listed  above  has  there  been 
any  indication  of  eddy  coalescence.  Papailiou  &  Lykoudis 
noted  that  the  breakdown  of  the  large  eddies  is  accomplished 
through  the  separation  of  the  eddies  into  smaller  rotating 
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parts.  These  parts  initially  remained  within  the  vortex 
but  eventually  separated,  forming  yet  smaller  lumps  of 
rotating  fluid.  Taneda  observed  that  this  breakdown  is 
often  followed  by  a  regeneration  of  the  two-dimensional 
structure  at  a  point  farther  downstream.  The  concept  of 
generation,  growth  and  decay  cycles  for  the  large  structures 
has  been  put  forward  by  Townsend.  His  recent  investigation 
provides  some  evidence  of  the  phenomena  and  he  has  illus¬ 
trated  the  degradation  to  correlation  measurements  that 
results  from  averaging  over  eddies  in  different  stages  of 
their  life  cycle. 

Bevilaqua  &  Lykoudis  and  Papailiou  &  Lykoudis  have 
observed  that  entrainment  in  plane  wakes  occurs  in  a  manner 
similar  to  that  observed  by  Dimotakis  &  Brown  in  the  mixing 
layer.  In  the  first  stage,  ambient  fluid  is  entwined  into 
the  large  rotating  eddies  by  turbulent  shearing.  The 
second  stage  of  complete  assimilation  of  the  recently 
entrained  fluid  is  accomplished  at  the  point  of  eddy  break¬ 
down  as  described  above. 

As  in  the  plane  mixing  layer,  there  is  some  question 
as  to  the  degree  of  two-dimensionality  of  the  large  eddies. 
Barsoum,  Kawall  &  Keffer  (1978)  have  made  measurements  of 
the  space-time  correlations  derived  from  the  intermittency 
function  for  a  plane  wake.  They  found  that  the  bulges  in 
the  turbulent-nonturbulent  interface  have  a  characteristic 
streamwise  dimension  approximately  double  the  spanwise 
dimension.  In  addition,  there  was  no  evidence  of  periodicity 
of  the  intermittency  function  with  respect  to  downstream 
distance.  Both  of  these  features  are  incompatible  with  the 
concept  of  an  ordered  vortex  street  structure. 

3.  Jet  Flows 

The  majority  of  the  reported  research  into  the  large 
scale  structures  of  jet  flow  has  dealt  with  the  initial 
region.  This  is  the  result  of  the  importance  of  these 


structures  to  noise  production  in  that  zone.  The  following 
investigations  are  representative  of  this  type:  Bradshaw 

et  al.  (1964),  Becker  &  Massaro  (1968),  Wooldridge  &  Wooten 
(1971),  Ko  &  Davies  (1971),  Crow  &  Champagne  (1971),  Lau, 
et  al.  (1972),  Fuchs  (1972),  Liu  (1974),  Merkine  &  Liu 
(1975),  Lau  &  Fisher  (1975),  Bruun  (1977),  Ff owe s  Williams 
&  Kempton  (1978) .  In  the  initial  region,  the  flow  field  is 
essentially  that  of  a  plane  mixing  layer  wrapped  around  the 
lip  of  the  nozzle.  Ring  vortices  are  shed  from  the  lip 
more  or  less  periodically  and  go  through  a  process  of  growth 
through  repeated  eddy  coalescence  as  they  are  convected 
downstream.  Due  to  the  inherent  instabilities  of  a  vortex 
ring,  the  structures  become  increasingly  three-dimensional 
so  that  by  the  end  of  the  potential  cone,  an  apparently 
fully  turbulent  flow  is  established.  Noise  production  in 
the  initial  region  appears  to  be  controlled  by  the  initially 
highly  coherent  ring  vortices  and  the  intense  first  few 
eddy  coalescence  events. 

Comparatively  fewer  experiments  have  been  undertaken 
in  the  similarity  region  of  turbulent  jets.  The  strongest 
indication  of  an  ordered  structure  is  the  well  documented 
apparent  "flapping"  behavior  of  turbulent  two-dimensional 
jets.  This  characteristic  is  easily  detected  by  measuring 
the  correlation  between  the  fluctuating  longitudinal 
velocities  from  opposite  sides  of  the  jet  centerline.  For 
a  zero  time  delay  the  result  is  a  negative  correlation 
coefficient  indicating  a  lateral  motion  of  the  instantaneous 
velocity  profile.  In  the  past,  this  behavior  has  been 
compared  to  the  flapping  of  a  flag  in  which  the  centerline 
of  the  mean  profile  is  periodically  displaced  about  its 
average  location.  It  was  thought  that  the  flapping  was 
essentially  a  traveling  wave  instability  rather  than  an 
indication  of  organized  large  scale  structure.  Flapping 
two-dimensional  jets  have  been  reported  by  Bradbury  (1965) , 


16 


Goldschmidt  &  Bradshaw  (1973),  Weir  &  Bradshaw  (1975), 

Gutmark  &  Wygnanski  (1976),  Everitt  &  Robins  (1978)  and 
Cervantes  (1978) . 

In  the  detailed  investigation  by  Cervantes,  the 
frequency  of  the  flapping  was  inferred  from  the  oscillating 
correlation  function  that  resulted  from  introducing  a  time 
delay  to  one  of  the  signals.  The  frequency  thus  determined 
was  found  to  be  independent  of  the  centerline  to  probe 
distances,  thereby  indicating  a  structure  which  spanned  the 
entire  flow.  A  similar  measurement  of  the  autocorrelation 
function  for  the  fluctuating  transverse  velocity  component 
at  the  centerline  led  to  the  same  flapping  frequency. 

Cervantes  also  found  that  the  variation  of  the  flapping 
frequency  with  respect  to  downstream  distance  is  compatible 
with  the  similarity  scaling  relationships  for  the  plane  jet. 
Consequently,  the  flapping  frequency,  when  nondimension- 
alized  by  the  local  mean  velocity  and  length  scales,  is 
equal  to  a  constant,  Sf  =  ffb/Um  *  0.1.  Approximate  agreement 
with  this  value  may  be  inferred  from  the  investigations  of 
Bradbury,  Everitt  &  Robins,  Goldschmidt  &  Bradshaw  and  Weir 
&  Bradshaw.  The  similarity  scaling  and  agreement  between 
the  various  investigators  on  the  dimensionless  flapping 
frequency,  Sf,  implies  a  Reynolds  number  independence  that 
was  verified  by  Cervantes. 

While  the  traveling  wave  phenomena  is  a  plausible 
explanation  for  the  flapping  jet  behavior,  an  alternate 
possibility  is  the  presence  of  large  scale  eddies  of  a 
von  KSrmSn  vortex  street  type.  Although  more  order  is 
implied  by  this  structure  than  traditionally  expected  for  a 
turbulent  plane  jet,  the  similarity  between  the  character¬ 
istics  of  the  jet,  wake  and  mixing  layer  make  it  a  plausible 
hypothesis  (see  Figure  2-5).  In  particular,  the  weaving 
flow  evident  in  the  flapping  jet  is  strongly  suggestive  of 
the  trailing  wake  vortex  street.  The  similarity  scaling 
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Figure  2-5  Organized  Structures  in  Two-dimensional 

Free  Shear  Flows 
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and  Reynolds  number  independence  of  the  jet  structural 
frequency  are  also  two  of  the  most  distinctive  features  of 
the  large  eddies  in  the  plane  mixing  layer. 

Recent  flow  visualization  experiments  by  Moallemi 
(1980)  strongly  suggest  a  vortex  street  type  pattern  of 
the  large  structures  in  a  plane  jet.  By  stretching  oil 
coated  smoke  wires  parallel  to  and  on  either  side  of  the 
jet  centerline  in  the  ambient  fluid,  the  silhouette  of 
apparently  two-dimensional  large  eddies  could  be  discerned 
as  the  smoke  was  carried  into  the  turbulent-nonturbulent 
interface.  Stretching  a  single  wire  normal  to  the  center- 
line  revealed  that  the  entrainment  velocity  vector  has  a 
negative  longitudinal  component  that  is  approximately  15% 
of  the  lateral  component. 

Mattingly  &  Criminale  (1971)  have  made  a  spatial 
stability  analysis  of  the  two-dimensional  plane  jet.  They 
found  that  the  most  unstable  and  rapidly  amplified  mode  led 
to  the  generation  of  a  vortex  street  type  of  structure.  It 
was  also  illustrated  that  this  structure  would  be  the  most 
effective  in  the  transfer  of  energy  from  the  mean  flow  to 
the  turbulent  field. 

In  summary,  a  simplified  hypothetical  model  of  the 
overall  structural  development  of  the  plane  jet  may  be 
described  in  terms  of  the  three  regions. 

1.  In  the  initial  region,  small  perturbations  of  the 
initially  laminar  boundary  layers  are  amplified  into  the 
characteristic  two-dimensional  vortices  of  a  turbulent 
mixing  layer.  The  characteristic  scales  in  this  region  are 
the  centerline  velocity,  UQ,  and  momentum  thickness,  Sj* 

The  mixing  layers  grow  until  at  a  downstream  distance  which 
scales  on  the  slotwidth,  they  contact  each  other  and  a 
fundamental  change  in  the  structure  takes  place. 

2.  Within  the  transition  region,  an  overall  pattern  of 
the  large  eddies  evolves  which  produces  instantaneous  and 


periodically  repeatable  antisymmetric  velocity  profiles. 

The  length  and  velocity  scales  within  this  region  include  the 
halfwidth  and  mean  centerline  velocity. 

3.  Further  downstream,  the  jet  has  adjusted  itself  to 
the  new  eddy  patterns  and  a  state  of  dynamic  equilibrium 
and  structural  similarity  is  established.  In  the  work  that 
follows,  the  similarity  region  will  be  the  focus  of 
attention. 


CHAPTER  III 


EXPERIMENTAL  TECHNIQUES 


A.  Flow  Field  Apparatus 

All  results  presented  herein  were  obtained  from 
measurements  in  the  two-dimensional  plane  jet  depicted  in 
Figure  3-1.  The  flow  was  generated  by  a  3/4  hp  blower 
discharging  to  a  56  x  51  cm  plenum  chamber.  A  gradual 
contraction  led  to  a  15.2  x  30.4  cm  channel  where  the  flow 
passed  through  straighteners  and  screens  before  the  final 
contraction  to  a  1.27  x  30.4  cm  slot  and  discharging  normal 
to  a  91.4  x  30.4  cm  wall.  The  developing  jet  was  bounded 
above  and  below  by  91.4  x  120.0  cm  parallel  walls.  Heating 
for  the  jet  was  accomplished  with  a  4  kW  mesh  wire  heating 
element  in  the  plenum  chamber. 

The  jet  was  operated  at  a  constant  Reynolds  number 
4  -1 

(Re,,  =  1.74  x  10  ,  =  21  msec  )  and  overheat  (A©  = 

12  °C)  above  the  ambient  temperature. 

Mean  flow  measurements  were  made  with  a  pitot  tube 
while  fluctuating  velocity  and  temperature  signals  were 
derived  from  the  output  of  a  TSI  1050  series  anemometer. 

The  flow  temperature  was  equal  to  the  ambient  while  velocity 
measurements  were  made  so  that  there  was  no  temperature 
contamination  of  the  signal.  An  overheat  of  1.8  in  the 
constant  temperature  mode  was  used  for  those  measurements. 
Fluctuating  temperature  measurements  were  made  using  a 
constant  current  probe  with  I  *  1.5  mA.  An  analysis  by 
Jenkins  (1974)  with  the  same  apparatus  revealed  that  the 
maximum  contamination  error  of  the  temperature  signal  due 


to  the  velocity  fluctuations  is  less  than  4%. 

B.  Turbulence  Detection  and  Intermittenc 


Measurement 

The  existence  of  a  sharply  defined  interface  sepa¬ 
rating  the  turbulent  and  nonturbulent  flow  near  the  free 
stream  boundary  of  a  turbulent  shear  flow  was  first 
reported  by  Corrsin  (1943)  for  an  axisymmetric  circular  jet. 
The  intermittent  or  burst-like  behavior  of  a  hot  wire  anemo¬ 
meter  signal  resulting  from  the  alternate  convection  of 
turbulent  and  nonturbulent  flow  past  the  probe  was 
quantified  by  Townsend  (1949) .  He  defined  the  intermittency 
function  I(r,t)  as  a  square  wave  which  takes  on  the  value  of 
unity  or  zero  corresponding  to  periods  of  turbulence 
(vortical  flow)  or  nonturbulence  (nonvortical  flow) ,  respec¬ 
tively,  at  a  fixed  point  in  the  flow  field.  The  intermittency 
Y (r) ,  may  be  found  by  taking  the  time  average  of  I(r,t)  and 
is  interpreted  as  the  fraction  of  time  that  the  flow  is 
turbulent  at  a  particular  position.  The  time  derivative  of 
I(r,t)  yields  a  train  of  alternating  positive  and  negative 
unit  impulses  that  correspond  to  the  unique  events  of  lead¬ 
ing  and  trailing  edge  interface  crossings.  The  frequency 
of  occurrence  of  the  positive  (or  negative)  impulses  is 
defined  as  the  interface  crossing  frequency,  f^(r). 

Townsend  also  developed  an  intermittency  detection 
instrument  that  became  the  predecessor  to  many  similar 
devices  that  have  followed.  The  primary  features  of  these 
circuits  are: 

1)  Definition  of  a  turbulence  indicator  signal,  S(r,t), 
that  is  derived  from  some  measurable  flow  property  which 
exhibits  a  strong  contrast  between  nonturbulent  low  level 
portions  of  the  signal  and  turbulent,  higher  level  segments. 

2)  Application  of  a  gate  or  threshold  level  to  S(r,t) 
which  defines  an  intermediate  square  wave,  G(r,t),  taking  on 
the  values  of  one  or  zero  corresponding  to  intervals  for 
which  S(r,t)  is  above  or  below  the  threshold  C. 
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3)  Application  of  a  hold  time,  t.  ,  to  G(r,t)  where  t. 

n  n 

is  the  prescribed  minimum  length  of  time  during  which 
— ►  t 

I(r,t)  is  allowed  to  make  two  changes  of  state.  The  hold 
time  provides  a  means  of  discriminating  between  the  minimum 
burst  (turbulent  segment  of  S(r,t))  and  void  (nonturbulent 
segment  of  S(r,t))  time  intervals  expected  from  the  flow 
and  the  apparent  shorter  ones  resulting  from  high  frequency 
turbulent  fluctuations  of  S(r,t)  and  scales  of  interface 
motion  below  the  range  of  interest. 

4)  The  analysis  technique,  typically  a  choice  between  an 
analog  circuit  or  computer  analysis  of  digitized  data 
records. 

The  choices  made  for  S(r,t),  C,  and  the  analysis 
technique  are  described  below.  A  representative  review  of 
the  literature  relating  to  turbulence  detection  and  inter- 
mittency  measurement  is  given  in  Appendix  A. 

1 •  The  Indicator  Signal.  S(r.t) 

The  indicator  signal  used  for  the  present  investiga¬ 
tion  is  the  instantaneous  temperature.  Turbulence  detection 
is  based  on  the  discrimination  between  slightly  heated 
turbulent  fluid  and  relatively  uniform  temperature  ambient 
fluid.  This  type  of  technique  has  been  utilized  by  several 
previous  investigators  (Barsoum,  et  al.  (1977),  Freymuth  & 
Uberoi  (1973),  Hedley  &  Keffer  (1974),  and  Kohan  (1969))  and 
a  systematic  evaluation  of  the  technique  is  given  by  LaRue 
(1974)  . 

The  primary  advantage  of  temperature  measurement  as  a 
criterion  for  S(r,t)  is  that  the  signal  provides  a  strong 
contrast  between  turbulent  bursts  and  nonturbulent  voids 
without  encountering  the  zero  crossing  problem.  This  is 
illustrated  in  Figure  3-2  which  is  a  photograph  of  an 
oscilloscope  trace  of  S(r,t).  Since  only  positive  deviations 
from  the  ambient  temperature  level  occur,  rectification  is 
unnecessary  and  the  number  of  extraneous  threshold  crossings 


3-2  Oscilloscope  Trace  of  S(r,t) 
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is  much  less  than  for  S(r,t)  based  on  velocity  fluctuations 
or  their  derivatives.  This  characteristic  allows  a  good 
degree  of  confidence  in  the  choice  of  C  and  since  the 
dependence  of  the  measured  intermittencies,  crossing 
frequencies  and  conditional  averages  on  these  parameters  is 
diminished. 

All  measurements  described  in  this  document  were  made 
with  a  jet  overheat  of  12°C  at  the  exit.  This  allowed  a 
centerline  mean  temperature  excess  of  3  to  8°C  over  the 
measurement  range  which  is  comparable  with  the  4.5°C  used  by 
Barsoum,  et  al.  (1978),  and  considerably  less  than  the  15° 
to  22°C  range  used  by  Jenkins  and  Goldschmidt  (1974) .  Mean 
velocity  and  temperature  measurements  were  taken  to  verify 
that  they  were  consistent  with  the  similarity  principles  for 
a  plane  jet.  Repetition  of  the  measurements  at  various 
spanwise  locations  did  not  reveal  any  appreciable  deviations 
of  the  mean  flow  from  two  dimensionality  over  the  region  of 
interest  due  to  buoyancy  effects. 

2 .  The  Threshold,  C 

The  role  of  the  threshold  in  the  turbulence  detection 
procedure  is  illustrated  in  Figure  3-3.  In  choosing  the 
optimum  C,  it  is  desirable  to  set  it  as  close  as  possible  to 
the  ambient  level  without  encroaching  into  the  noise  band 
surrounding  the  ambient  signal.  The  noise  referred  to  here 
is  believed  to  be  the  result  of  small  amplitude  temperature 
fluctuations  in  the  room  air.  If  the  threshold  is  chosen 
below  the  maximum  level  of  the  noise  fluctuations,  false 
indications  of  interface  crossings  will  result. 

There  are  two  consequences  of  being  unable  to  use  the 
ideal  threshold,  that  is,  C  infinitesimally  above  the 
ambient  level.  The  first  of  these  is  that  the  intermittency 
function  will  be  systematically  underestimated  due  to  the 
finite  slope  of  the  leading  and  trailing  edge  temperature 
profiles  as  shown  in  Figure  3-3.  The  second  is  that  the 
probability  of  mistaking  a  short  period,  low  temperature  dip 
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within  a  turbulent  burst  as  a  true  separation  between 
bursts  increases  with  increasing  C.  This  is  also  due  to 
the  finite  rise  and  fall  time  of  the  temperature  level 
since  the  dips  are  of  increasingly  longer  duration  for 
higher  thresholds.  Since  it  is  the  function  of  the  hold 
time  to  compensate  for  the  short  period  dips  below  C,  a 
higher  threshold  increases  the  uncertainty  in  the  choice  of 
This  problem  has  also  been  discussed  by  LaRue  (1974). 

The  choice  of  threshold  for  the  present  investigation 
was  made  from  subjective  inspections  of  S(r,t)  for  various 
intermittencies.  From  these,  the  width  of  the  noise  band 
around  the  ambient  level  could  be  estimated  and  the  thresh¬ 
old  set  accordingly.  Visual  observations  suggested  that  the 
amplitude  of  the  noise  band  did  not  vary  significantly  with 
longitudinal  or  transverse  position.  Since  the  primary 
consideration  in  choosing  C  was  the  width  of  the  noise  band, 
no  advantage  would  be  gained  from  varying  C  at  different 
locations.  To  have  varied  C  would  have  in  fact  introduced 
additional  uncertainty  into  the  measurement  of  y  and  f^  due 
to  their  systematic  dependence  on  C  as  described  in  the 
previous  paragraph  .  The  dependence  of  and  f^  on  the 
threshold  is  illustrated  in  Oler  &  Goldschmidt (1978) . 

3-  The  Hold  Time.  _Tlt 

The  function  of  the  hold  time  is  to  eliminate  short 
period  "apparent"  voids  within  bursts  or  inversely,  bursts 
during  voids  as  illustrated  in  Figure  3-3.  As  pointed  out 
by  Kibens  and  Kovasznay  (1970) ,  an  investigation  of  the 
intermittency  characteristics  of  a  particular  flow  field 
should  be  predicated  on  a  definition  of  the  degree  of  detail 
desired  concerning  the  interface  geometry.  For  the  present 
investigation,  the  large  scale  eddies  are  of  principal 
interest. 

The  nominal  value  of  at  x/D  =  30  was  determined 
(using  the  preselected  threshold)  by  arbitrarily  choosing 
the  value  which  would  give  agreement  between  the  measured 
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crossing  frequency  and  the  frequency  as  determined  from 
visual  inspection  of  strip  chart  recordings  of  S(r,t) (at 
the  position  of  maximum  crossing  frequency) .  The  value 
determined  in  this  way  at  x/D  =  30  was  x^  «  4  msecs.  It 
was  held  constant  for  all  measurements  in  the  transverse 
plane  at  that  x/D. 

For  x/D  /  30,  the  hold  time  was  varied  in  direct 

proportion  with  a  global  time  scale  derived  from  the 

variation  of  the  centerline  velocity  and  velocity  halfwidth. 

For  a  self-preserving  plane  jet,  the  halfwidth,  b,  is 

expected  to  vary  linearly  witi  x/D^and  the  centerline 

velocity,  U^,  should  vary  as  ^  .  This  behavior  is 

demonstrated  in  Figures  2-3  and  2-4.  The  variation  of  the 

p  ”13 /2 

ratio  b/Um  is  expected  to  go  as  ^  and  may  be  inter¬ 
preted  as  a  global  time  scale  for  the  plane  jet.  The  hold 
time  variation  was  proportional  to  this  scale  with  the 
constant  of  proportionality  determined  at  x/D  =  30.  This 
resulted  in  .7  msec  <  xh  <  11.5  msec  for  10  <  x/D  <.  60  which 
is  consistent  with  the  3.6  -  7.5  msec  range  utilized  by 
Jenkins  and  Goldschmidt  (1974),  for  25  <.  x/D  <  55  in  the 

same  flow  field.  Their  choice  of  x,  was  based  on  a 

h  ■* 

characteristic  frequency  derived  from  S(r,t)  as  suggested  by 

Kibens  and  Kovasznay  (1970) .  Ajagu  and  Goldschmidt  (1976) 

used  a  constant  x.  =  10  msec  for  25  <  x/D  <.  45  for  their 

h 

measurements  in  a  plane  jet.  In  addition,  Kohan  (1969) 
chose  3.4  msec  <  x^  <.  4.1  msec  in  a  two-dimensional  wall  jet 
and  Kibens  and  Kovasznay  (1970),  used  x^  =  4  msec  for  a 
turbulent  boundary  layer. 

The  sensitivity  of  the  interface  crossing  frequency  to 
xh  is  illustrated  in  Figure  3-4.  Variation  of  the  hold  time 
produces  a  large  change  in  the  measured  frequencies  and  a 
distortion  of  the  normalized  distribution  of  f^/f^.  These 
effects  are  the  result  of  varying  the  degree  of  interface 
detail  included  in  the  formulation  of  the  intermittency 
function. 
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The  analogous  case,  in  which  is  held  constant  and 
the  time  scales  of  the  interface  motion  change  for 
increasing  x  is  illustrated  in  Figure  3-5.  Distributions 
of  fY^fym  measure<*  at  x/D  =  30,  40,  and  50  exhibit  approxi¬ 
mately  the  same  discrepancies  as  in  Figure  3-4.  Clearly, 
if  a  meaningful  comparison  of  interface  motion  is  to  be  made 
over  a  wide  range  of  streamwise  positions,  the  hold  time 
must  be  varied  appropriately.  Failure  to  do  so  results  in 
an  inconsistent  discrimination  between  the  various  scales 
of  interface  motion  in  the  formulation  of  the  intermittency 
function. 

4 .  The  Analysis  Technique 

Although  the  exact  manner  in  which  S,  C  and  have 
been  incorporated  varies  widely  in  the  literature,  the  basic 
analysis  techniques  may  be  classified  as  either  analog  or 
digital.  Analog  intermittency  detection  instruments 
typically  offer  low  set-up  and  operating  costs.  In 
addition,  there  is  a  certain  amount  of  convenience  that 
results  from  having  a  continuous,  real  time  intermittency 
signal  available  for  analysis.  These  features  are  obtained 
at  the  expense  of  some  precision  in  the  application  of  C  and 
Th  to  S(r,t)  to  obtain  I(r,t).  Digital  intermittency  analy¬ 
sis  techniques  offer  a  high  degree  of  precision  that  is 
limited  only  by  the  resolution  and  speed  of  the  analog-to- 
digital  converter.  This  precision  is,  however,  obtained 
with  much  higher  equipment  purchase  and  operating  costs  than 
for  a  similar  analog  instrument.  In  addition,  digital 
analysis  techniques  are  typically  less  conveniently  accom¬ 
plished  in  the  laboratory  environment  than  analog  techniques. 
Any  improvement  in  convenience  in  digital  system  usage  is 
usually  associated  with  much  higher  costs. 

Two  methods  of  digital  analysis  have  been  utilized  in 
the  present  research.  In  the  early  stages  of  the  program, 
analog  recordings  of  the  indicator  signal  were  made  for 
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playback  and  analysis  using  a  General  Automation  mini¬ 
computer.  The  details  of  the  procedure  and  the  fortran 
programs  involved  are  described  in  detail  by  Oler  & 
Goldschmidt  (1978)  and  in  Appendix  B. 

The  final  analysis  procedure  which  was  utilized  for 
the  majority  of  the  intermittency  measurements  utilizes  an 
original  hybrid  technique.  The  application  of  the  threshold 
to  the  indicator  signal  was  accomplished  with  an  analog 
comparator  circuit.  Application  of  the  hold  time  to  this 
signal  was  then  performed  by  an  online  microcomputer  whose 
output  was  the  completely  conditioned,  real-time  intermit¬ 
tency  signal.  A  description  of  the  analog  circuitry,  micro¬ 
computer  and  machine  language  code  is  given  by  Oler  & 
Goldschmidt  (1979)  and  in  Appendix  C. 

With  the  intermittency  signal  available  as  a  continuous 
real-time  analog  signal,  it  was  necessary  to  measure  the 
average  value  (intermittency)  and  the  pulse  rate  (interface 
crossing  frequency) .  The  former  was  accomplished  by  using 
the  intermittency  signal  to  control  the  gate  of  a  digital 
counter  whose  input  was  a  10  kHz  sine  wave.  Allowing  the 
counter  to  run  over  an  interval  of  100  sec  with  gate 
intervals  specified  by  the  intermittency  function  resulted 
in  a  final  count  which  was  some  fraction  of  1  x  10^.  This 
fraction  is  equivalent  to  the  intermittency  fraction.  The 
interface  crossing  frequency  was  measured  by  inputing  the 
intermittency  function  directly  to  the  counter  and 
averaging  over  100  sec. 

C .  Correlation  Measurements 

Correlation  measurements  made  during  the  research 
program  consisted  of  two  types:  1)  two  point  fluctuating 
longitudinal  velocity  cross-correlations  and  2)  two  point 
intermittency  function  cross-correlations. 

1.  Velocity  Cross-correlations 

All  fluctuating  velocity  correlations  were  measured 
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utilizing  a  Thermo- Systems,  Inc.  analog  correlator  and 
true  rms  voltmeters.  Due  to  the  low  normalized  correlation 
coefficients  which  were  involved,  it  was  necessary  to  use 
the  100  sec  time  constant  settings  and  allow  four  time 
constants  for  settling. 

2 .  Intermittent  Cross-correlations 

For  normalized  correlations  it  is  necessary  to  have 
signals  with  zero  mean  values.  To  correlate  intermittency 
signals,  it  was  not  possible  to  AC  couple  or  high  pass  filter 
the  signals  prior  to  input  to  the  analog  correlator  due  to 
the  distortion  of  the  square  wave  patterns  and  lack  of 
accuracy  that  would  result.  An  original  alternate  method 
was  utilized  which  took  advantage  of  the  logic  level  pulse 
train  character  of  the  intermittency  signals. 

The  normalized  correlation  coefficient  is  defined  as: 
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_  lim  .  T/2 

(  )  =  t-*-®  =  f  (  )  dt 

-T/2 


Y1  =  Il(t)  =  I(rl't) 


Y2  =  I2(t)  =  I(r2,t) 


Utilizing  the  fact  that  =  3*1  may  be  expanded  as 

I.  (t)  I-  (t  +  T)  -  Yn  Y2 

R^ Ar ,  T )  =  =4 - - -  —  -1-  (3.2) 


{(I1(t)2  -  Y12][I2(t  +  T)2  -  y22I } 1/2 
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Eq.  3.2  may  be  simplified  further  utilizing  the  fact  that  I(t) 
only  takes  on  the  values  (0,1).  Therefore, 

Ht)2  =  I(t)  =  Y 
so  that  Eq.  3.2  becomes 


Rj (Ar, r) 


I1tt)I2(t  +  t)  -  y1Y2 

Uy^T^Uy^YpT)1^ 


(3.3) 


The  intermittencies  in  Eq.  3.3  were  determined  using  the  gated 
oscillator  technique  described  previously.  Noting  that  the 
product  of  two  logic  level  signals  is  the  equivalent  of  the 
two  signals  logically  ANDed  together,  the  same  technique  was 
utilized  to  measure  the  time  average  of  I1(t)I2(t  +  t)  after 
performing  the  logical  AND  operation  electronically. 

An  additional  complication  in  measuring  the  cross¬ 
correlation  between  two  intermittency  signals  resulted  from 
the  on-line  microcomputer  being  capable  of  applying  the  hold 
time  to  only  one  signal  at  a  time.  This  problem  was  remedied 
by  applying  the  hold  time  to  one  signal  and  recording  the 
output  simultaneously  with  the  other  signal.  The  hold  time 
could  be  applied  to  the  second  signal  on  playback  resulting 
in  the  availability  of  both  intermittency  signals  for  the 
zero  time  delay  correlations.  This  procedure,  utilizing  the 
recorder,  has  the  additional  feature  of  allowing  the  time 
averages  1^ (t)  ,  I2 (t  +  x)  and  1^ (t) I2 (t  +  t)  to  be  taken 
over  simultaneous  data  records. 

Time  delays  were  introduced  using  a  second  on-line 
microcomputer.  It  was  programmed  to  sample  the  fully  formu¬ 
lated  intermittency  functions  at  5  kHz,  store  the  result  in  a 
data  buffer  and  simultaneously  output  the  saved  status  of  the 


intermittency  signals  with  a  prescribed  delay.  The  machine 
language  code  required  for  this  function  is  given  in 
Appendix  C. 
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CHAPTER  IV 
EXPERIMENTAL  RESULTS 


The  primary  goal  of  the  research  program,  as  noted  in 
Chapter  I,  was  to  investigate  the  coherent  structures 
of  a  fully  developed  plane  jet.  Early  in  the  research,  the 
presence  of  a  vortex  street  type  structure  was  suspected  and 
hence  hypothesized  for  the  remainder  of  the  work.  The  task 
then  became  to  test  the  compatibility  of  the  vortex  street 
model  with  respect  to  measurable  jet  properties. 

A  vortex  street  type  structure  would  have  to  exhibit: 

1)  Periodicity 

2)  Organization  over  large  distances 

3)  Some  degree  of  two-dimensionality 

4)  Compatibility  with  similarity  considerations. 

The  requirement  for  two-dimensionality  should  not  be 

interpreted  too  rigorously.  Conceivably,  a  vortex  structure 
could  consist  of  a  complex,  yet  organized  and  periodic 
structure  such  as  that  observed  by  Perry  &  Lim  (1978)  and 
Pao  &  Kao  (1977) .  Although  their  investigations  were  in 
axisymmetric  fields,  the  possibility  of  a  slanted  or  looped 
structure  in  a  two-dimensional  field  (especially  one  confined 
by  horizontal  planes)  should  not  be  discarded. 

Compatibility  with  similarity  considerations  is  a 
requirement  if  the  structure  is  to  play  a  consistent  role  in 
the  transport  dynamics  of  a  self-preserving  flow.  However, 
if  the  structure  were  simply  propagated  or  convected  through 
by  the  mean  field,  similarity  would  not  be  necessary  for  its 
existence.  As  it  turns  out,  every  indication  has  been  that 
the  primary  source  of  energy  transfer  from  the  mean  flow  to 
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the  fine  scale  turbulence  is  the  large  scale  structures. 

From  that  point  of  view,  structural  similarity  would  be 
required  by  the  concept  of  local  dynamic  equilibrium  in  a 
developing  and  self-preserving  flow.  In  Chapter  V,  the 
inducing  nature  of  the  vortex-like  structure  will  be 
discussed,  showing  that  it  is  not  simply  a  wave-type  motion 
impressed  on  the  mean  flow. 

The  major  objective  of  the  experimental  program  has 
been  to  evaluate  the  extent  of  compliance  of  the  real  jet 
flow  with  the  vortex  street  characteristics  described  above. 
The  properties  measured  were  the  intermittency,  interface 
crossing  frequency  and  space- time  correlations.  These  and 
inferences  from  them  are  now  presented. 

A.  Intermittency  and  Interface  Crossing 
Frequency 

A  vortex-like  structure  should  dominate  the  character¬ 
istics  of  the  flow  field  and  dictate  the  instantaneous 
position  of  the  turbulent-nonturbulent  interface.  If  the 
geometry  of  the  interface  is  controlled  by  the  underlying 
turbulent  structure,  then  properties  of  that  structure  may 
be  inferred  from  intermittency  measurements. 

The  lateral  distributions  of  the  intermittency, 
y,  and  normalized  interface  crossing  frequency,  f^/f^m,  are 
given  in  Figures  4-1  and  4-2.  The  lateral  coordinate  used 
for  these  plots  has  been  centered  at  the  position  of  maximum 
crossing  frequency  (ym/b  =  1.6  +  .05)  and  scaled  by  the 
local  velocity  halfwidth,  b.  For  this  jet  (see  Figure  2-4), 

b/D  =  .098  (x/D  +  .619)  (4.1) 

when  x/D  >  10 

The  plots  exhibit  similarity  in  that  the  curves  from 
various  longitudinal  positions  collapse  on  the  y/b 
coordinate. 


The  mean  interface  position  is  given  by  the  point  of 
maximum  crossing  frequency.  At  that  location,  as  noted  from 
Figure  4-1,  the  corresponding  intermittency  is  approximately 
0.6,  not  0.5  as  would  be  demanded  by  a  symmetric  distribution 
function.  The  same  feature  (y  >  0.5  at  f  =  f  )  has  been 
noted  by  Gutmark  &  Wygnanski  (1976),  Mourn,  et  al.  (1979)  and 
Everitt  &  Robins  (1978).  The  explanation  is  that  the  inter¬ 
face  is  more  sharply  folded  in  its  troughs  than  on  its 
crests.  In  the  present  investigation,  this  feature  indicates 
that  the  interface  maintains  a  qualitatively  similar  shape 
at  all  streamwise  positions  since  the  y  =  0.6  value  is 
maintained. 

The  fact  that  the  curves  collapse  on  a  y/b  plot  also 
implies  similarity  for  the  underlying  lateral  turbulent 
length  scales  and  geometries.  Similarity  scaling  or  simply 
structural  similarity  is  critical  to  proving  the  validity 
of  the  vortex  street  model  since  the  large  eddies  are 
expected  to  play  a  consistent  role  in  the  downstream 
development  of  the  jet.  Additional  indications  of  structural 
similarity  will  be  noted  throughout  the  rest  of  this  chapter 
and  the  concept  will  provide  the  key  to  formulation  of  the 
interpretive  model  given  in  the  next  chapter. 

The  downstream  variation  of  £^m  is  given  in  Figure  4-3 
along  with  the  apparent  "flapping"  frequency  of  Cervantes 
(1978) .  Cervantes  measured  the  frequencies  with  the  cross¬ 
correlation  function  of  longitudinal  velocity  fluctuations 
from  opposite  sides  of  the  jet.  f^  corresponds  to  the 
frequency  inferred  from  the  time  delay  for  the  first  peak  of 
the  oscillating  correlation  function  while  f2  represents 
that  derived  from  the  average  of  the  peak  to  peak  delay 
intervals.  Rather  than  "flapping"  frequency,  the  term 
structural  passage  frequency  (or  simply  structural  frequency) 
will  be  used  throughout  the  remainder  of  this  report. 


Strong  support  for  the  assumption  of  relatedness 
between  the  interface  geometry  and  underlying  structure  is 
given  by  Figure  4-3.  Although  both  experiments  were  made 
at  the  R.  W.  Herrick  Laboratories,  they  were  completely 
independent  of  each  other  and  were  made  on  totally 
separate  jet  apparatus. 

From  similarity  arguments  for  the  plane  jet,  global 

length  scales  are  expected  to  vary  in  proportion  with  x/D 

-1/2 

while  velocity  scales  should  go  as  (x/D)  .  Therefore,  a 


-3/2 


global  frequency  scale  should  vary  as  (x/D) ‘ .  This 

behavior  is  verified  by  the  crossing  frequency  data  of 

Figure  4-3.  If  a  dimensionless  frequency  is  formed  from 

the  structural  frequency,  f  ,  velocity  halfwidth  and  mean 

centerline  velocity,  a  flow  constant,  S  ,  results.  Values 

s 

of  this  constant  for  the  present  research  and  reported 

measurements  of  "flapping"  frequency  are  given  in  Table  4.1. 

The  relatively  close  agreement  suggests  a  universal 

instability  mechanism  similar  to  that  analyzed  by  Mattingly 

&  Criminale  (1971) (see  Section  II. B. 4). 

From  the  streamwise  variation  of  f  ,  it  may  be  shown 

ym 

that  the  net  number  of  large  scale  structures  per  unit 
length  decreases  as  the  flow  develops.  If  u  (x/D)  is  the 

5 

convective  velocity  of  the  structures,  then  the  average 
wavelength  of  the  pattern  is  given  by 


i  (x/D)  =  u  (x/D) /f  (x/D) 

S  SS 


(4.2) 


Since  the  pattern  represented  by  Figure  2-5  consists  of  two 
vortices  per  wavelength,  the  density  of  vortices  in  a  cross¬ 
stream  section  of  unit  thickness  is 


n (x/D)  =  2/1 _(x/D) 
s 


Combining  Eqs.  4.2  and  4.3  yields 


Table  4-1.  Summary  of  Reported  Structural  Frequencies 


fs  b 


Um 

x/D 

rd 

Investigator 

0.1  -  0.244" 

44-71 

- 

Bradbury  (1965) 

0.06* 

45 

2.6  x  104 

Goldschmidt  & 
Bradshaw  (1973) 

0.12* 

12 

105 

Weir  &  Bradshaw 
(1975) 

0.08  -  0.14+ 

20-100 

7.9  x  103 
-  1.5  x  104 

Cervantes  (1976) 

o.i  -  o.is44" 

100 

3  x  104 

Everitt  &  Robins 
(1978) 

.09* 

10-60 

4 

1.7  x  10 

present 

*  Determined  from  average  peak  to  peak  intervals  of  the 
correlation  function  between  velocity  fluctuations 


+  Determined  by  interval  to  the  first  peak  of  the  correlation 
function  between  velocity  fluctuations 

^  Determined  by  peak  in  cross-spectral  density  distribution 
of  velocity  fluctuations 

Determined  from  interface  crossing  frequencies 
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n (x/D)  =  2fs(x/D)/us(x/D)  .  (4.4) 

Treating  the  structural  density  as  a  transferable  flow 
quantity,  a  measure  of  the  net  local  "production"  rate  of 
the  vortices  is 


qn(x/D)  =  Us(x/D)  (4.5) 

The  variation  of  u  ,  f  and  n  may  be  described  by  the 

s  s 

similarity  scaling  relationships  given  in  Eqs.  4.6 


ug(x/D) 

__ 


=  Cy (x/D) 


-1/2 


Us  =  V 


-1/2 


(4.6a) 


f  (x/D) D 

-  =  Cp  (x/D) 


“3/2  •  F  =  C  X‘3/2 
'  s 


C  C 

n (x/D) D  (X/D)'1  ;  N  =  2  X_1 

‘■'U 


Substituting  Eqs.  4.6  into  Eq.  4.5  yields 


(4.6b) 

(4.6c) 


O 


0 

=  q 


n 


D2/U 


o 


=  U  dN/dX 
s 

=  -2  CpX”5/2  (4.7) 


From  Figure  4-3,  C_  =  2.09  and  Eq.  4.7  becomes 

r 
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Qn  =  -4.18  X_5/2  (4.8) 

Although  a  net  reduction  in  the  number  of  large  structures 
is  indicated  by  the  minus  sign  in  Eg.  4.8,  no  conclusions 
may  be  drawn  concerning  the  mechanism  which  accounts  for  that 
loss.  Either  eddy  coalescence  or  a  net  unbalance  of  the 
natural  generation  and  decay  of  the  large  eddies  are  viable 
possibilities. 

B.  Spatial  Characteristics  of  the  Large  Structures 

A  vortex  street  structure  in  a  fully  developed  jet  flow 
would  have  to  exhibit  specific  spatial  relationships  for  the 
instantaneous  velocity  field  and  turbulent- nonturbulent 
interface  geometry.  In  particular,  the  counter  rotating 
vortices  positioned  on  opposite  sides  of  the  centerline 
would  produce  antisymmetric  transverse  distributions  of 
longitudinal  velocity  and  intermittency .  Additionally, 
there  would  be  a  periodically  repeated  pattern  of  the 
velocity  field  and  interface  geometry  with  respect  to  down¬ 
stream  distance. 

The  existence  of  the  spatial  relationships  described 
above  may  be  tested  experimentally  by  measurements  of  space- 
time  correlations  of  velocity  fluctuations  and  the  inter¬ 
mittency  function.  The  correlation  coefficients  of  the 
longitudinal  velocity  fluctuations  and  intermittency 

functions  at  points  A  and  B  with  separations  5  ,  6  and  6 

x  y  z 

and  relative  time  delay,  t,  are  given  by  Eqs.  4.9  and  4.10. 


vt,vvv  = 


UA*t*UB*t+X* 


UA  V 


(4.9) 
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RI(Tf6x'6y'6z) 


(iA(t)-y(iB(wl‘ifB| 

>(V''a)2(ib-V2i1/2 


(4.10) 


1 .  Velocity  Correlation  Results 

The  instantaneous  asymmetry  of  the  longitudinal 
velocity  distribution  is  evidenced  by  the  distributions  of 
Ru(0;°»*5y*0)  given  in  Figure  4-4.  The  most  obvious  feature 
of  the  distributions  are  the  strong  negative  correlations 
between  the  velocity  fluctuations  on  opposite  sides  of  the 
centerline  and  extending  well  into  the  intermittent  regions. 
Similar  results  have  been  found  by  Goldschmidt  &  Bradshaw 
(1973),  Gutmark  &  Wygnanski  (1976),  Cervantes  (1978)  and 
Everitt  &  Robins  (1978). 

The  correlation  measurements  which  result  from  both 
longitudinal  and  lateral  separation  of  the  velocity  probes 
are  presented  in  Figure  4-5  in  the  form  of  an  isocorrelation 
contour  map.  The  map  provides  a  graphic  illustration  of  the 
instantaneously  antisymmetric  velocity  field  and  periodically 
repeatable  structural  pattern.  The  actual  distance  over 
which  the  measurable  flow  organization  exists  (x/D  =  20  to 
45)  is  quite  remarkable  and  strongly  supports  the  possibil¬ 
ity  of  a  vortex  street  eddy  pattern. 

2 .  Intermittency  Correlation  Results 

The  correlation,  R ^  (0 ; 0 , 5^ , 0 ) ,  was  measured  for  the 
intermittency  signals  at  points  symmetrically  placed  on 
either  side  of  the  centerline  such  that  Y.  =  yq  =  0.5.  A 
vortex  street  structure  should  force  an  antisymmetric  dis¬ 
placement  of  the  interfaces  such  that  the  resulting  corre¬ 
lation  is  negative.  The  measurements  revealed,  however, 
a  zero  correlation,  indicating  that  the  motions  of  the 
interfaces  were  independent  of  one  another.  The  same 
conclusion  was  found  by  Wygnanski  &  Gutmark  (1971)  and 
Mourn,  et  al.  (1979) . 


m 


tk 
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The  periodicity  of  a  vortex  street-type  structure 
should  produce  periodic  interface  geometries.  An  indication 
of  this  should  be  present  in  autocorrelations  of  the  inter- 
mittency  functions.  Results  of  measurements  of  this  type 
are  given  in  Figure  4-6  (x/D  =  30,  y  =  0.5).  For 
comparison,  similar  experimental  results  by  Mourn,  et  al. 
are  included  (x/D  =  10,  y  =  0.5).  For  a  meaningful  compari¬ 
son  of  the  experiments,  the  time  delay,  t,  has  been  made 
dimensionless  by  the  local  length  and  velocity  scales,  b 

and  U  ,  respectively, 
m 

Although  there  is  no  indication  of  periodicity  in  the 
present  investigation,  the  results  of  Mourn,  et  al.,  do 
contain  a  weak  negative  lobe.  They  ascribe  the  negative 
correlations  to  statistical  instabilities  that  result  from 
applying  Fourier  transforms  to  finite  length  data  records. 

If  instead,  the  time  interval  for  the  negative  peak  were 
interpreted  as  a  half-period  of  the  interface  motion,  a 
dimensionless  structural  frequency  of  Sg  =  0.06  would  be 
indicated.  This  interpretation  of  the  Mourn,  et  al.  data 
compares  favorably  with  the  values  of  Table  4.1,  supporting 
further  the  hypothesis  of  a  vortex  street  structure. 
Unfortunately,  the  present  data  do  not  exhibit  the 
expected  negative  lobe.  This  may  be  attributed  to  either 
the  different  techniques  used  in  the  determination  of 
I(r,t)  and  Rj.(t),  to  the  difference  in  measurement  positions 
(x/D  =  30  vs.  x/D  =  10)  or  to  violation  of  the  hypothesis. 

C .  Two-Dimensionality 

The  vortex  street  structure  of  a  fully  developed  plane 
jet  flow  should  exhibit  a  certain  degree  of  two-dimensional¬ 
ity,  subject  to  the  qualifications  given  at  the  beginning  of 
this  chapter.  This  two-dimensionality  should  be  reflected 
in  the  internal  turbulent  velocity  field  as  well  as  in  the 
interface  geometry. 

Quantification  of  the  degree  of  two-dimensionality  is 
accomplished  experimentally  by  comparing  characteristic 


longitudinal  and  vertical  length  scales.  These  length 
scales  are  derived  from  the  correlations  of  velocity  and 
intermittency  (see  Eqs.  4.9  and  4.10)  and  are  defined  as: 


Ax  =  0/O°  V0;6x,0,6z)d6x  (4.11) 

Az  =  0/O°  Rg{0;6x'0'5z)d6z  (4.12) 


with  g  =  u  or  I. 

1 .  Length  Scales  of  the  Intermittent  Region 

Figure  4-7  illustrates  correlations  of  the  inter¬ 
mittency  function  at  x/D  =  30  and  y  =  0.5.  Similar  measure¬ 
ments  were  taken  for  10  <  x/D  <60.  It  should  be  noted  that 
the  space  correlations,  Rj-tOjfi  ,0,0),  like  the  autocorre¬ 
lations,  Rj(t;0,0,0)/  of  the  previous  section,  do  not 
exhibit  a  negative  lobe.*  (A  negative  lobe  would  have 
further  confirmed  the  hypothesized  periodicity  of  the  inter¬ 
face  geometry) . 

The  integral  length  scales  were  derived  from 
Rt(0;6  , 0,0)  and  R  (0;0,0,6  )  curves  such  as  Figure  4-7 

A  X  Z 

using  Eqs.  4.11  and  4.12. 

The  resulting  streamwise  distributions  of  A1  and  A1 

X  z 

are  given  in  Figure  4-8.  Once  again,  the  self-preserving 
variation  of  all  intermittency  length  scales  is  evident. 

From  Figure  4-8,  it  is  apparent  that  the  longitudinal  and 
vertical  length  scales  are  approximately  equal. 


*  Auto  and  space  correlations  in  the  mean  flow  direction  may 
be  related  by  invoking  Taylor's  hypothesis,  i.e., 

R(o;<5  ,o,o)  =  R(-u  t, -o,o,o)  . 


b 
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AX/D  =  .047  (x/D  +  2.5)  (4.13) 

or 

AX/b  -  .480  (4.14) 


Mourn,  et  al.  have  also  reported  the  near  equivalency 

of  A1  and  A1.  At  x/D  =  10,  they  report  A*  -  A*  =  .47b. 

X  Z  X  z 

2.  Length  Scales  of  the  Turbulent  Velocity  Field 

Determination  of  the  integral  length  scales  at  x/D  = 

30  and  y/b  =  0.0,  0.25,  0.5  was  accomplished  utilizing 
correlations  of  the  longitudinal  velocity  fluctuations.  In 
this  case,  the  integration  with  respect  to  the  separation 
distances  was  limited  to  the  intervals  over  which  the  corre¬ 
lations  were  positive. 

A  typical  result  is  plotted  in  Figure  4-9.  The  average 
length  scales  determined  from  the  correlation  curves  are 


A^/b  =  .65  +  .03  (4.15) 

A"/b  =  .28  +  .02  (4.16) 

Z 

The  indication  here  is  that 

AU/AU  «  2  ,  (4.17) 

X  z 

which  is  surprising  considering  the  measurements  in  the 
intermittent  region.  However,  the  results  are  consistent 
with  those  reported  by  Everitt  &  Robins  (1978),  i.e.. 


A^/b  =  .47  +  .02  (4.18) 

AU/b  =  .23  +  .02  (4.19) 

z  — 


indicating 


A^/A*  *  2  . 

X  Z 


(4.20) 


Another  unexpected  feature  of  Figure  4-10  is  the 
negative  lobe  in  the  R(0;0,  0,  6z)  correlation.  There  are 
three  possible  explanations: 

1)  A  combination  of  streamwise  tilting  of  the  large 

vortices  and  the  streamwise  periodicity  of  the 

velocity  field  (see  Figure  4-5)  results  in  the 

negative  R  (0;0,0,6  )  correlations, 
u  z 

2)  The  vortices  undulate  periodically  along  their  length 
and  the  negative  lobe  in  the  correlation  is  an 
indication  of  the  first  wavelength. 

3)  There  is  a  lack  of  coherence  of  the  vortices  in  the 
vertical  direction  so  that  the  correlation  is 
dominated  by  the  near  isotropic  small  scales.  The 
negative  lobe  is  then  the  result  of  continuity  require¬ 
ments  in  isotropic  turbulence  correlations. 

Based  on  the  comparison  of  integral  scales  alone,  the 
degree  of  two-dimensionality  characteristic  of  the  large  struc¬ 
tures  is  quite  limited.  However,  two  points  must  be  noted: 

1)  The  flow  visualization  experiments  of  Moallemi  (1980) 
have  indicated  an  apparently  two-dimensional  structure. 

2)  A  similar  contradiction  between  visualization  and 
inferences  from  integral  scales  concerning  two- 
dimensionality  is  also  present  in  the  plane  mixing 
layer.  (See  Wygnanski  &  Fiedler  (1970) ,  Chandrsuda, 
et  al.  (1979)  ,  Pui  &  Gartshore  (1979)  and  Wygnanski, 
et  al.  (1979) )  . 

D.  Convective  Velocities 
In  this  section,  four  different  experimental 
approaches  to  the  determination  of  the  convective  velocity 
of  the  vortex-like  structure  are  presented. 

The  classicial  technique  for  the  determination  of 
convective  velocities  is  to  determine  the  space-time  corre¬ 
lation  function  for  velocity  probes  separated  in  the  mean 
flow  direction.  By  noting  the  separation  distance  and  the 
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time  delay  which  yields  the  peak  correlation  coefficient, 
the  convective  velocity  may  be  determined  from  their  ratio. 
Measurements  of  this  type  have  been  made  by  Ott  (1972)  and 
the  results  are  tabulated  in  Table  4.2.  A  variation  on  this 
approach  is  to  band-pass  filter  the  velocity  signals  prior 
to  forming  the  correlation  function.  Young  (1973)  made 
measurements  utilizing  this  technique  and  those  are  also 
given  in  Table  4.2. 

Between  the  measurements  of  Young  and  Ott,  those  by 
Young  based  on  the  lowest  frequencies  should  be  most  closely 
tied  to  the  velocity  of  the  underlying  large  scale  struc¬ 
ture.  However,  the  difficulty  with  the  results  of  both 
experiments  is  that  the  convective  velocities  decrease  mono- 
tonically  away  from  the  centerline  (but  not  as  a  fixed  ratio 
of  the  local  mean  velocity) .  There  is  no  clear  choice  for  a 
single  velocity  that  could  be  associated  with  the  underlying 
pattern. 

Cervantes  (1978)  measured  a  convective  velocity  of 
the  "flapping"  motion  which,  based  on  the  comparisons  of 
Section  IV. A,  should  correspond  to  the  convective  velocity  of 
the  vortex-like  large  structures.  In  addition  to  the 
symmetric  spacing  of  the  probes  on  opposite  sides  of  the 
centerline,  Cervantes  also  introduced  a  streamwise  separation. 
By  observing  the  shift  of  the  negative  peak  correlation,  the 
convective  velocity  was  estimated  in  a  manner  analogous  to 
that  of  Young  and  Ott. 

The  experiment  was  repeated  for  two  values  of  the  lat¬ 
eral  spacing  and  the  results  are  given  in  Table  4.2.  Again, 
the  convective  velocities  are  noted  to  be  nonconstant  and 
their  relationship  to  a  unique  structural  velocity  is  not 
clear. 

In  the  presently  reported  research  program,  a  convec¬ 
tive  velocity  was  determined  two  ways. 

First,  it  was  obtained  from  space-time  correlations  of 
the  intermittency  function.  Examples  of  these  measurements 
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Table  4-2.  Summary  of  Structural  Related  Convective 
Velocities 


WWW 

uc^um 

y/b 

Investigator 

1.0 

.0 

+Ott  (1972) 

.95 

.5 

.82 

.75 

.68 

1.0 

.57 

1.5 

.425 

0 

Young  (1973) 

.39 

.5 

.39 

.75 

.38 

1.0 

.30 

1.5 

.75 

.65 

t 

Cervantes 

.61 

1.0 

.27 

1.6 

++present 

** 

.50 

- 

present 

+  Measured  from  broad  band  space-time  correlations 

*  Measured  from  lowest  frequency  narrow  band  correlations 

^  Measured  by  extension  of  flapping  frequency  correlations 

Measured  from  space-time  correlations  of  the  intermittency 
function 

**  Measured  utilizing  flow  visualization  technique 

***  The  symbol  ' u  '  is  used  to  represent  all  of  'convective' 
velocities  described  IV. D 
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for  various  probe  separations  are  given  in  Figures  4.10  and 
4.11.  The  results  are  given  in  Figure  4-12  and  were 
tabulated  in  Table  4.2.  The  indicated  convective  velocities 
of  the  interface  bulges  are  uCI/um  ~  0.27  +  .02. 

Secondly,  in  collaboration  with  Moallemi  (1980),  the 
convective  velocity  of  the  structural  pattern  was  obtained 
from  the  instantaneous  positions  of  the  vortices  using  flow 
visualization  and  consideration  of  the  self-preserving 
nature  of  plane  jets. 

At  a  given  longitudinal  position,  the  structural 
frequency,  pattern  wavelength  and  structural  velocity  are 
related  by  Eq.  4.2,  or 

u  =  £  *  f  (4.21) 

s  s  s 

Since  the  jet  is  a  self-preserving  flow,  Eq.  4.21  may  be 
expressed  in  terms  of  the  similarity  variables: 


U  =  L  *  F 
s  s  s 


where  ,  following  Eqs.  4.6 


Fs  "  CFX 


-3/2 


Ls  =  CLX 


Us  "  CUX 


-1/2 


Substitution  of  Eqs.  4.23  into  Eq.  4.22  yields 


CU  -  CL  *  CF 


(4.22) 


(4.23a) 

(4.23b) 

(4.23c) 


(4.24) 


If  CL  and  Cp  are  determined  experimentally,  then  the  Cy  will 
be  given  by  Eq.  4.24  and  the  structural  velocity  by  Eq.  4.23c. 
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From  Section  IV. A  (Table  4.1),  the  structural 
frequency  constant  was  defined  as 

fsb 

Sc  =  -rf-  =  0.1.  (4.25) 

s  Um 
m 

For  the  plane  turbulent  jet  (see  Figures  2-3  and  2-4) 


b/D  =  0.1  (x/D) 

U  /U  =2.36  (x/D)~1/2 
m  o 

Combining  Eqs.  4.26  and  4.27  with  Eq.  4.25  yields 

f  D 

Fs  = 

_  D  Um 
"  s  b  UQ 


*  2.36 


X~  3/2 


Therefore,  from  Eq.  4.23a, 


(4.26) 

(4.27) 


(4.28) 


Cp  =  2.36 


(4.29) 


In  the  following  chapter,  it  will  be  shown  that  the 
longitudinal  positions  of  the  vortices  in  a  vortex  street 
which  obeys  the  similarity  laws  for  a  plane  jet,  are 
uniquely  defined.  In  particular,  if  a  vortex  is  instanta- 

4*  Vi 

neously  at  and  the  position  of  the  n  downstream  vortex  is 

Xd,  then  (see  Eq.  5.27) 
n 


C 


L 


2  ln(Xn/XQ) 
n 


(4.30) 


As  a  consequence  of  a  Eq.  4.30,  if  the  positions  of  n+l 
vortices  are  indicated  from  a  photograph,  then  n  estimates 
of  CT  may  be  determined. 

L 
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The  above  procedure  was  followed  using  the  flow 

visualization  results  of  Moallemi.  It  was  found  that  due  to 

the  finite  downstream  length  required  for  self-preservation 

to  be  established,  C_  is  a  function  of  X  that  approaches  its 

L 

similarity  value  asymptotically.  Therefore,  the  value  of 
C  defined  by  (Xn,X  )  using  Eq.  4.30  was  assigned  to  an 
intermediate  position,  X,  representing  the  equi-distant 
point  in  time  between  XQ  and  X^  and  given  by  (see  Chapt.  V) 

~CLn 

X  «  XQe  r  (4.31) 

The  results  of  this  experiment  are  presented  in 
Figure  4-13.  The  asymptotic  value  found  is 

C  =  0.5  (4.32) 

Li 

Substitution  of  this  value  for  CL  and  Cp  from  Eq.  4.29 
into  Eq.  4.24  yields 

Cy  =  1.18  (4.33) 


From  Eq.  4.23c, 


U  =  1.18  X_1/2  (4.34) 

s 

The  ratio  of  the  structural  velocity  and  mean  centerline 
velocity  comes  from  Eqs.  4.27  and  4.34,  giving 

u  /U  =0.5  (4.35) 

s  m 

This  result  was  also  tabulated  in  Table  4.2. 

A  review  of  Table  4.2  shows  a  substantial  difference 

between  the  broadband  U  and  the  other  values.  This  is 

c 
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understandable  since  it  includes  the  high  frequency 
components  not  characteristic  of  the  large  scales.  The 
extrapolated  values  of  Young,  those  now  measured  from 
intermittency  correlations  and  those  inferred  from  photo¬ 
graphs  are  not  too  different  from  one  another.  Those  from 
Cervantes  "flapping"  measurements  are  also  in  agreement 
although  somewhat  larger  in  magnitude.  All  of  these  are 
based  on  different  global  characteristics,  each  of  which  is 
attributed  to  the  vortex-like  structure.  The  differences 
noted  could  possibly  be  explained  by  the  inherent  differences 
in  the  measurement  processes. 

The  highest  confidence  is  placed  upon  the  visualization 
data  as  a  representation  of  the  actual  velocity  of  the 
vortex  cores.  The  corresponding  value  of  Us/Um  =  0*5  will 
be  used  in  the  model  derivation  in  the  following  chapter. 

E.  Summary 

From  the  experimental  tests  of  the  validity  of  a 
vortex  street  model  for  the  large  scale  structure  of 
turbulent  plane  jets,  the  following  conclusions  are  found: 

1)  The  global  variation  of  all  length,  velocity  and 
frequency  scales  are  compatible  with  similarity 
requirements. 

2)  The  large  vortex-like  motion  is  periodic. 

3)  The  coherent  vortex-like  motion  extends  across  the 
entire  width  of  the  flow  and  for  considerable  up  and 
downstream  distances. 

4)  The  structural  organization  inferred  from  the 
measurements  has  the  instantaneously  antisymmetric 
characteristic  expected  of  a  vortex  street  structure. 

5)  Only  a  limited  degree  of  two-dimensionality  may  be 
inferred  from  the  measurements. 

The  structural  pattern  is  convected  at  a  velocity  near 
one-half  the  mean  centerline  velocity. 


6) 
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CHAPTER  V 

A  NUMERICAL  SIMULATION 


A.  Introduction 

As  an  alternate  means  of  testing  the  compatibility  of 
a  vortex  street  model  with  the  known  characteristics  of 
plane  jet  flows,  a  numerical  simulation  has  been  performed. 
The  term  simulation  is  chosen  to  distinguish  this  approach 
from  one  which  makes  predictions  of  turbulent  properties 
based  upon  the  solution  of  a  set  of  modeling  equations. 
Instead,  the  simulation  provides  an  investigative  tool  with 
which  a  kinematic  representation  of  the  structure  may  be 
evaluated.  The  validity  of  the  hypothesized  structure  is 
judged  by  the  extent  of  agreement  between  the  calculated 
flow  properties  and  experimental  and  theoretically  based 
expectations. 

The  approach  taken  to  a  representation  of  the  hypoth¬ 
esized  vortex  street  structure  is  that  of  a  linear  super¬ 
position  of  Rankine  vortices.*  The  velocity  at  any  point 
in  the  flow  field  is  then  determined  by  the  additive 
effects  of  all  vortices  making  up  the  pattern.  Although 
a  reasonable  approximation  might  have  been  to  superpose 
the  vortex  street  structure  onto  an  otherwise  homogeneous 
mean  flow,  that  approach  was  not  utilized  in  the  presently 
reported  research.  It  was  assumed  instead  that  the 
vortices,  themselves,  are  totally  responsible  for  the  flow 


*  The  results  of  identical  calculations  using  a  slightly 
different  model  for  the  individual  vortices  are  given 
in  Appendix  D. 
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field.  The  validity  of  this  assumption  will  be  discussed 
in  the  final  section  of  this  chapter. 

The  key  consideration  in  the  formulation  of  the  model 
is  the  strict  application  of  similarity  scaling  requirements. 
Experimental  results  given  in  the  preceding  chapter 
repeatedly  indicated  the  characteristic  of  structural 
similarity  in  plane  jet  flows.  That  property  is  utilized 
to  prescribe  the  streamwise  variation  of  all  the  geometric 
and  dynamic  parameters  which  define  the  hypothesized 
structure. 

Another  feature  of  the  simulated  structure  which  is 
notable  for  its  absence,  is  the  effect  of  the  small  scale 
turbulence.  In  the  current  research  program,  the  large 
scale  structures  have  been  of  principal  interest  and 
consequently,  only  those  effects  have  been  modeled.  The 
question  of  how  the  large  structures  interact  with  the 
small  scales  is  a  difficult  problem  in  itself.  Undoubtedly, 
the  small  scale  turbulence  is  critically  important  to  some 
of  the  calculated  flow  properties  and  the  consequences  of 
their  neglect  is  evident.  However,  by  isolating  the 
effects  of  the  large  scales,  it  is  clearly  indicated  which 
of  the  other  flow  properties  are  dominated  by  those  scales. 

Similar  kinematic  representations  of  the  vortex 
street  in  a  two-dimensional  near  wake  have  been  reported  by 
Hooker  (1936)  and  Schaefer  &  Eskinazi  (1959).  These 
investigations  utilized  potential  vortices,  Rankine 
vortices  or  a  combination  of  both  to  simulate  the  stable 
formation  region  of  the  wake.  The  calculations  were  not 
extended  into  the  fully  developed  turbulent  field. 

Takaki  &  Kovasznay  (1978)  have  utilized  a  physical 
representation  of  the  large  vortices  in  a  temporally 
developing  mixing  layer  to  predict  the  eddy  coalescence 
rate.  Their  primary  interest  was  in  the  prediction  of  the 
statistical  distribution  of  vortex  spacings  and  as  such, 
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did  not  attempt  to  calculate  the  velocities  induced  by  the 
vortices. 

To  the  author's  best  knowledge,  the  structural  model 
formulation  described  herein  is  the  first  to  fully  account 
for  the  consequences  of  global  turbulent  flow  field 
development  on  the  determination  of  local  mean  properties. 
As  will  be  shown,  this  is  the  key  to  understanding  the 
effects  produced  by  the  coherent,  large  scale  structures. 

B.  Formulation  of  the  Model 

The  basic  elements  that  make  up  the  numerical 
simulation  of  the  vortex  street  model  are:  1)  specifica¬ 
tion  of  the  velocity  induced  by  an  isolated  vortex  in  terms 
of  its  characteristic  length  and  velocity  scaling  para¬ 
meters,  2)  description  of  the  basic  vortex  street  geometry 
and  specification  of  the  variation  of  its  defining  para¬ 
meters  with  respect  to  downstream  position,  3)  development 
of  a  rationale  for  the  determination  of  the  longitudinal 
coordinates  of  each  vortex  making  up  the  pattern,  and  4) 
calculation  of  the  flow  properties  from  the  net  effect  of 
the  entire  vortex  street. 

1 .  The  Rankine  Vortex 

The  vortex  street  structure  is  modeled  with  an  array 
of  Rankine  vortices.  These  vortices  have  the  classical 
circumferential  velocity  field  of  a  diffusing  line  vortex 
(see  Lamb  (1945))  given  by 

v6  2tTt  11  ■  exp(“Ivt)1  (5,1) 

and  illustrated  in  Figure  5-1.  Evaluation  of  d  \^/dr  =  0 
reveals  that  the  maximum  velocity  occurs  at 

r  =  core  radius 
c 

/T726 

"/  4  vt  .  (5.2) 


Substitution  of  Eq.  5.2  into  Eq.  5.1  yields 


|1  -  exp  (-1 . 26-^-*-)  ]  . 

r_ 


(5.3) 


From  Eq.  5.3  and  Fig.  5-1,  it  is  noted  that  the 
velocity  field  surrounding  a  Rankine  vortex  is  defined  by 
the  parameters  r,  the  circulation  strength  and  rc»  the  core 
radius.  For  r  >  the  induced  velocity  is  essentially 

that  of  a  potential  vortex  and  for  r  <  r  ,  there  is 

C 

approximately  solid  body  rotation. 

If  a  Rankine  vortex  is  centered  at  (x^,y^)  in  the  two- 
dimensional  coordinate  system  shown  in  Figure  5-2,  the 
velocity  components  at  (x,y)  are  given  by 


^(x.y)  = 


-rtnt  r.  2 

1  ^(1  -  exp t- 1.2 6^)  ]} 

ci 


27rri 


r.C,  r.  2 

vi(x,y)  =  2i7^2  t1  ~  exp[-1.26(~~)  ]} 


where 


(5.4) 


(5.5) 


r 


2 


i 


=  x 


xi 


■  y  - 


The  parameters  which  determine  the  induced  velocities  at 
(x,y)  and  whose  streamwise  variation  must  be  specified 
from  similarity  scaling  relationships  are 


Lateral  Velocities 
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>1  "  *i<xi> 


(5.6) 


°i  ‘  rci(xi> 


(5.7) 


Fi  =  W 


(5.8) 


2 .  Similarity  Scaling  of  the  Vortex  Street 

The  hypothesized  vortex  street  representation  of  a 

plane  jet  is  illustrated  in  Figure  5-3.  From  the  figure, 

it  is  noted  that  the  defining  parameters  of  the  pattern  at 

any  longitudinal  station  are  the  vortex  center  offset 

distance,  h(x),  the  core  radius,  r  (x) ,  and  the  circulation 

c 

strength,  T  (x)  . 

Also  indicated  in  the  figure  is  the  indexing  system 
that  will  be  used  throughout  this  chapter.  One  of  the 
vortices  with  a  negative  lateral  location  is  arbitrarily 
chosen  as  the  reference  position  and  labeled  (xQ,yo)  •  Its 
properties  are  then  given  by 


yo  =  -h(xQ) 


(5.9) 


V  =  rc(xo> 
o 


(5.10) 


Fo  =  "r(V  ' 


(5.11) 


The  properties  of  the  i  up  or  downstream  vortex  are  given 
by 


yA  =  (-Di+1  h(Xi) 


(5.12) 


rc.  =  W 

r.  =  (-i)i+1  r(x.) 


(5.13) 


(5.14) 


where 


d 

x  .  =  x . 
i  1 


u 

or  x^ 


Since  the  plane  jet  is  a  self-preserving  flow  and 
structural  similarity  is  expected  for  the  vortex  street, 
the  variation  of  h,  rc  and  r  are  defined  by  the  global 
length  and  velocity  scales.  In  particular. 


h (x/D) 
D 


rc(x/D) 

D 


r(x/D) 

U  D 
o 


=  CH  *  (g)  ;  H (X)  =  CHX 

"  Cr*  »  RC(X)  “  crX 

=  cr*  (£)  /  ?  r(x)  =  crx1/2 


(5.15) 

(5.16) 

(5.17) 


Substituting  Eqs.  5.15  -  5.17  into  Eqs.  5.12  -  5.14  yields 


it 

•H 

>* 

(-1)1+1  H(Xi) 

(5.18) 

R 

Ci 

w 

(5.19) 

Fi  = 

(-D1+1  r(x.) 

(5.20) 

where 

X.  = 

l 

u  d 

Xi  or  Xi  . 

3.  Specification  of  the  Longitudinal  Vortex  Coordinates 
The  key  feature  remaining  undefined  is  the  instanta¬ 
neous  longitudinal  coordinate,  xV  or  X^  ,  of  each  vortex. 
As  the  vortices  move  downstream,  the  pattern  wavelength 
increases,  the  convective  velocity  decreases  and  the  number 
of  vortices  per  unit  length  diminishes.  All  of  these 
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features  must  be  accounted  for  in  the  determination  of 
vortex  positions  utilizing  a  rationale  that  is  compatible 
with  similarity  scaling  requirements.  This  has  been 
accomplished  in  the  presently  reported  formulation  by 
utilizing  the  concept  of  vortex  density  introduced  in 
Sections  IV. A  and  IV. B.  Vortex  density  is  useful  in  that 
it  allows  the  description  of  a  delta  function  (the  vortex 
locations)  in  terms  of  a  continuous  function. 

At  a  longitudinal  station,  X,  the  vortex  pattern  wave¬ 
length  (see  Eq.  4.22)  is  given  by 

L  =  U  /F  (5.21) 

s  s  s 


Ls  =  CLX 


(5.22) 


It  should  be  noted  that  the  wavelength  Ls(X)  is  a  point 
property  since  length  scales  change  continuously  up  and 
downstream.  For  an  observer  fixed  at  X,  it  is  the  wave¬ 
length  that  he/she  would  calculate  from  knowledge  of  the 
pattern  velocity  and  period  at  that  point. 

A  similar  point  property  is  the  local  vortex  density 
(see  Eq.  4.3)  defined  by 


(5.23) 


The  constant  is  a  consequence  of  there  being  two  vortices 

in  each  wavelength  of  the  pattern. 

Consider  a  vortex  instantaneously  located  at  (Xo,Yq) 

which  is  arbitrarily  chosen  as  the  reference  vortex  as 

described  in  the  previous  section  and  depicted  in  Figure 

5-3.  The  number  of  vortices  between  X^  and  some  downstream 

o 


25 
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position  Xa  may  be  found  by  integrating  the  vortex  density 
with  respect  to  X.  This  quantity  is  defined  as  the  inte- 
grated  vortex  count,  M(XQ,X  ),  and  is  given  by 


M(X  ,XU)  =  v/  N  dX 
°  o 


-  x;  cHc  dx  * 

O  L 


(5.24) 


Similarly,  the  integrated  vortex  count  between  XQ  and  an 
upstream  position  Xu  is  given  by 


M  ( XQ  /  X  )  = 


dX  . 


(5.25) 


Given  the  reference  vortex  position,  X  ,  the  position 
of  the  next  downstream  vortex,  x“,  is  such  that 
M(Xq,X^)  =  1.  Utilizing  Eq.  5.24  yields 


M<x£,  xj)  =  1 


-  xf 

o 


ll  2 


therefore 


CL  X1 

T  “  ln  (5T> 


X1  =  xoexplTil 


(5.26) 


The  position  of  the  n1"  downstream  vortex  is  found  in  an 
analogous  manner: 


U 
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M(X  ,X*J) 
o  n 


n 


d  nW 

X°  =  xoexp[-^i]  (5.27) 

The  same  procedure  is  used  to  determine  the  position  of  the 
first  upstream  vortex  from  XQ, 

u  ”CL 

xj  =  Xoexp[— .  (5.28) 

A.U 

The  position  of  the  n  upstream  vortex  is  given  by 

u  ”n^'i 

x”  =  x^xpl-yi)  .  (5.29) 

To  summarize,  at  a  given  instant,  the  arbitrary  selec¬ 

tion  of  any  reference  vortex  (X  ,Y  ) (with  the  restriction 

o  o 

that  Yq  =-H (XQ) )  is  sufficient  to  completely  determine  the 
longitudinal  positions  of  all  other  up  and  downstream 
vortices  using  Eqs.  5.27  and  5.29,  respectively.  The 
corresponding  lateral  coordinates,  core  radii  and  circula¬ 
tion  strengths  may  then  be  found  from  Eqs.  5.18  -  5.20. 

4.  Calculation  of  Temporally  Averaged  Flow  Properties 

The  net  velocity  at  a  point  (X,Y)  resulting  from  the 
linear  superposition  of  Rankine  vortices  making  up  a 
vortex  street  is  found  by  a  simple  summation  of  the  effects 
of  the  individual  vortices.  Given  that  the  reference 
vortex  is  instantaneously  at  (X  ,Y  ) ,  the  resulting  instan¬ 
taneous  velocity  components  are  given  by  (see  Eqs.  5.4  and 
5.5) 


>u4fav* 
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where 


& 


=  Velocities  induced  by  the  reference  vortex 


at  <X0,Yo> 


=  Velocities  induced  by  the  upstream  vortex 


at  (xj,  xj) 


=  Velocities  induced  by  the  downstream  vortex 

at  (X*,X?) 


The  time  average  of  the  fluctuating  velocity  compo¬ 
nents,  U (X , Y , t)  and  V (x, Y, t) ,  may  be  found  by  integrating 

with  respect  to  time  over  one  period  of  the  pattern  as  it 

* 

passes  the  transverse  plane  at  X.  The  time  dependence  of 
U  and  V  will  then  be  determined  by  the  instantaneous 
positions  of  the  vortices  in  the  street.  When  performing 
this  operation,  provisions  must  be  made  for  the  fact  that 
equal  time  increments  will  not  result  in  equal  convection 
distances  of  the  vortices.  This  is  a  consequence  of  the 
streamwise  variation  of  length  and  velocity  scales. 

The  problem  is  again  resolved  by  utilizing  the  concept 
of  a  continuous  vortex  density.  Consider  the  convection  of 
vortices  past  the  measuring  station  X.  Over  the  time 
interval,  0  to  t,  the  number  of  vortices  which  pass  X  may 
be  defined  as  the  integrated  vortex  flux,  A(X,t).  This 
*  t  represents  a  dimensionless  time  ;  *t  =  tU^/D  • 
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quantity  is  given  by 


A(X,t)  =  0/fc  Ug  N  dt  (5.32) 

Noting  that  U  and  N  are  functions  of  X  only  and  using 

S 

Eq.  5.23,  Eq.  5.32  becomes 


A  (X ,  t) 


(5.33) 


Over  the  interval  of  one  period  of  the  motion  past  X, 
A (X, t)  will  increase  linearly  by  two,  indicating  the 
passage  of  two  vortices.  This  suggests  utilizing  the 
integrated  vortex  flux  as  a  means  of  determining  the  time 
dependent  positions  of  the  vortices  making  up  the  vortex 
street. 

The  shift  of  the  reference  coordinate,  XQ,  downstream 
during  the  integration  period  is  found  by  relating  the 
integrated  vortex  flux  and  integrated  vortex  count.  In 
particular,  if  for  XQ  =  X  at  t  =  0,  then  from  Eq.  5.27, 


A(X, t) 


2 

CLX 


dX 


or 

A(X,t)C. 

Xo  =  X  exp  [ - 5 - (5.34) 

With  the  instantaneous  position  of  the  reference  vortex 
given  by  Eq.  5.34,  the  longitudinal  coordinates  of  all 
other  up  and  downstream  coordinates  are  given  by  Eqs.  5.29 
and  5.27,  respectively.  This  is  sufficient  to  determine 
the  instantaneous  velocity  components  at  (X,Y;t)  so  that 
the  temporal  averages  can  be  found  in  a  conventional 


manner 
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C.  Details  of  the  Numerics 
The  vortex  street  model  formulated  in  the  previous 
section  is  based  upon  the  variation  of  four  vortex  related 
similarity  variables: 


L 

-  CLX 

(5.35) 

r 

=  CpX1'2 

(5.36) 

H 

-  CHX 

(5.37) 

Rc 

=  CRX 

(5.38) 

As  a  result,  there  are  four  proportionality  constants  that 
must  be  determined;  C L,  Cp,  CR  and  CR. 

The  choices  for  the  constants  have  been  based  upon 
experimental  estimates  or  optimization  of  the  calculated 
results.  The  length  coefficient,  CL,  was  determined  from 
the  flow  visualization  experiments  (see  Section  IV. D)  and 
set  equal  to  0.5.  Specification  of  the  circulation 
coefficient,  Cp,  was  accomplished  by  choosing  the  value 
which  gave  a  match  between  calculated  and  experimentally 
determined  mean  centerline  velocities  at  X  *  x/D  =  45  (i.e., 

U  =  U  /U  =  2.36*  (45) -;L/2 ) .  The  result  was  Cr  =  2.597. 

'  m  m  o  i 

The  values  chosen  for  the  offset  and  radius 
i  coefficients  optimize  the  comparison  between  experimental 

and  calculated  results.  Moallemi  (1980)  has  estimated 
j  CH  =  .04  from  his  visualization  experiments  while  CR  =  .05 

V  was  chosen  for  the  calculations.  C_  was  set  equal  to  .13 

:  ^  K 

or  ,26Cl.  The  primary  consideration  in  the  choice  of  CR 
and  C_  was  to  match  the  experimentally  determined 
spreading  rate,  i.e.,  b/D  =  0.1 (x/D). 

Table  V-l  provides  a  tabulation  of  the  coefficients 
which  have  been  utilized,  criteria  for  their  choice  and  a 


1  < 

i  * 
r 
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comparison  with  their  visually  estimated  counterparts. 

Of  the  four,  the  choices  for  CH  and  CR  are  the  most 
arbitrary.  Illustrations  of  the  sensitivity  of  the 
calculated  results  to  these  parameters  are  given  in 
Appendix  D. 

Additional  parameters  which  were  involved  in  the 
numerical  simulation  are  the  time  step  and  the  vortex  street 
truncation  length.  For  all  calculations,  the  time  interval 
corresponding  to  one  wavelength  of  travel  by  the  pattern  was 
divided  into  20  equal  time  steps.  Except  for  experiments 
concerning  the  truncation  length  specifically,  the  effects 
of  20  vortices  both  up  and  downstream  of  the  reference 
vortex  were  included  in  the  velocity  calculations.  For  a 
reference  vortex  at  x/D  =  45,  this  corresponded  to  a  stream- 
wise  range  of  0.64  <  x/D  <  5200  for  the  41  vortex  positions. 

D.  Numerical  Calculation  Results 
1.  Mean  Longitudinal  Velocities 

The  mean  longitudinal  velocity  at  (X,Y)  has  been 
calculated  as 


U(X, Y)  =  -  /T  U(X, Y,t)  dt  .  (5.39) 

T  o 

The  resulting  transverse  distributions  of  U(X,Y/B)/Um(X) 
are  illustrated  in  Figure  5-4  for  X  =  15,  30,  45,  60  and 
75.  Also  included  in  the  figure  are  analogous  experimental 
distributions  and  the  analytical  predictions  of  Gorter 
(1942)  and  Reichardt  (  1943). 

The  transverse  scaling  length,  or  halfwidth  (see 
Section  II. A. 3),  B,  for  Figure  5-4  as  well  as  all  the 
following  illustrations  of  lateral  distributions  is 
defined  as 


b/D  =  0  *  (x/D) 
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or 


B  =  0.1  *  X  .  (5.40) 

The  coefficient,  B,  is  defined  as  the  spreading  rate  and 
was  determined  experimentally  as  shown  in  Figure  2-4. 

The  collapse  of  the  calculated  profiles,  from  the 
various  transverse  planes  onto  the  reduced  coordinate,  Y/ B, 
supports  the  manner  in  which  the  self-preserving  character¬ 
istic  of  real  jet  flows  was  incorporated  into  the  formula¬ 
tion  of  the  model.  Transverse  distributions  of  all  other 
calculated  flow  properties  exhibit  the  same  degree  of 
collapse  and  consequently,  only  the  representative 
distribution  from  X  =  45  is  given  for  those  (tabulated 
results  will  be  given  in  Oler  &  Goldschmidt  (1980)). 

From  Figure  5-4,  it  is  noted  that  the  vortex  street 
model  accurately  represents  the  mean  longitudinal  velocity 
distribution  and  spreading  rate  &  determined  from  experi¬ 
mental  measurements.  To  further  verify  the  capability  of 
the  model  to  duplicate  similarity  scaling,  the  variation 
of  the  mean  centerline  velocity,  Um(X) ,  with  respect  to  X 
is  given  in  Figure  5-5.  From  Figure  2-3,  it  was  shown 
that 


Um(x/D)/U0  =  2.36  *  (x/D) -1/2 

U  (X)  =  C„  *  x"1/2  (5.41) 

m  u 

with  which  the  model  agrees.  It  should  be  recalled  that 

the  similarity  coefficient,  Cp,  was  chosen  such  that 

U  (45)  =  2.36  *  (45)~1/2.  The  coefficient  was  held 
m 

constant  for  all  other  calculations  and  the  mean  center- 
line  velocities  illustrated  in  Figure  5-5  are  the  results. 


ft 
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As  noted  in  Section  II. B. 3,  Moallemi  (1980)  has 
observed  that  the  velocity  vector  of  the  potential  flow 
beyond  the  jet  boundary  actually  has  a  negative  longitu¬ 
dinal  component.  Experiments  with  the  numerical  simulation 
have  demonstrated  that  this  phenomenon  may  be  a  consequence 
of  the  finite  axial  length  over  which  the  mean  velocity 
gradients  are  strong  enough  to  maintain  a  coherent  vortex¬ 
like  large  scale  structure. 

The  original  approach  taken  for  the  numerical 
simulation  was  to  extend  the  vortex  street  far  enough  down¬ 
stream  such  that  the  furthest  vortices  had  a  negligible 
effect  on  the  calculation  of  local  flow  properties.  The 
implicit  assumption  was  that  self-preservation  would  be 
maintained  for  0  <  X  <  00  .  While  theoretically  appealing, 
that  assumption  is  only  approximately  true  in  real  flows. 

Figure  5-6  illustrates  the  effect  of  a  finite  length 
of  the  vortex  street  and  Figure  5-7  magnifies  the  effect  on 
the  potential  outer  flow.*  It  is  apparent  that  the 
longitudinal  component  of  the  velocity  vector  approaches 
zero  from  negative  values  in  the  limit  as  the  streamwise 
length  of  the  vortex  street  goes  to  infinity. 

The  surprising  implication  of  these  results  is  that 
there  may  be  local  deviations  from  perfect  similarity  due 
to  the  influence  (or  rather  the  lack  of  it)  of  coherent 
structures  much  further  downstream. 

2.  Mean  Lateral  Velocities 

The  mean  lateral  velocities  were  calculated  numeri¬ 
cally  by 

*  The  numerical  simulation  does  not  provide  any  indication 
of  interface  position.  However,  the  maximum  extent  of 
the  intermittent  region  may  be  inferred  from  experi¬ 
mentally  based  expectations  of  its  relationship  to  the 
mean  velocity  profile. 
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V (X, Y,  t)  »  -  /T  V(X, Y,t)  dt  .  (5.42) 

T  ° 

Figure  5-8  illustrates  the  resulting  transverse 
distributions,  V(X,Y/B)/Um(X)  vs.  Y/B.  Again,  the  self¬ 
preserving  variation  of  length  and  velocity  scales 
calculated  by  the  model  are  evident. 

The  distribution  of  mean  lateral  velocity  may  also  be 
calculated  from  knowledge  of  the  mean  velocity  profile  and 
conservation  of  mass.  Continuity  requires 


(X, Y)  .  3V (X, Y)  _  A 

lx -  *  15 - 0 


V (X, Y)  =“/Y  —  dY 


From  similarity  for  the  plane  jet, 


(5.43) 


» <*•?/*>.  =  f  (y/b) 


vx> 


(5.44) 


V1X-Y/B)  -  g  (Y/B)  . 


U„(X) 

m 


(5.45) 


Substituting  Eqs.  5.44  and  5.45  into  Eq.  5.43  yields 


um(x)g(n)  =  B(x)o/"  If(n)um(x)j  dn 


(^)  =  -  ^ 
1  U 


B (X)  dUm(X)  ,  dBI 


f  («])d7 


dB  (X) 


where  r\  =  Y/B. 


7  £{V 


(5.46) 


92 


This  procedure  has  been  followed  by  Gutmark  &  Wygnanski 
(1976)  and  the  results  are  given  in  Figure  5-8  for  compari¬ 
son  with  the  distribution  calculated  by  the  vortex  street 
model.  From  the  agreement  of  the  two  analyses  evidenced  by 
Figure  5-8,  it  is  verified  that  the  simulated  flow 
satisfies  conservation  of  mass  requirements. 

Extension  of  the  integral  limit  of  Eq.  5.46  yields  the 
lateral  velocity  at  infinity  or  the  entrainment  velocity, 

V  (X) ,  i.e. 

6 


g(«<>) 


B  (X) 

Vx> 


dU  (X) 
m 

dX 


dB (X)  \ 

dX  0  J 


f  (<j)d»| 


(5.47) 


Townsend  (1976)  has  shown  that  g(°°)  *  .054  which  is  the 
asymptotic  value  indicated  in  Figure  5-8. 

3.  Reynolds  Stress  Distributions 

Local  values  of  the  Reynolds  stress  were  calculated 
by 

UV  (X ,  Y)  =  }■  /T  {(U(X,Y,t)  -  U(X,Y)  ]  [V(X,Y/t)  -V(X,Y)J)  dt 

m  O 

(5.48) 

_  —  2 

The  resulting  distributions  of  UV(X,Y)/Um(X)  are  given  in 
Figure  5-9. 

Also  given  in  the  figure  are  experimentally  determined 
distributions  of  the  Reynolds  stress  by  Jenkins  (1974)  and 
Gutmark  &  Wygnanski  (1976).  By  plotting  on  the  reduced 
coordinate,  Y/B,  the  distributions  exhibit  reasonable 
agreement. 

The  Reynolds  stress  distribution  may  also  be 
calculated  from  the  mean  velocity  profile,  continuity 
equation,  and  streamwise  component  of  the  Navier-Stokes 
equations.  With  the  appropriate  boundary  layer 
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approximations  for  plane  jet  flows  (see  Townsend  (1956)), 
the  Reynolds  averaged  X-component  of  the  Navier-Stokes 
equations  becomes 


U  (X,  Y) 


3U (X, Y) 

ax 


+  V(X,Y) 


au(xtY) 

3Y 


3UV  (X,  Y) 
3Y 


(5.49) 


Utilizing  the  continuity  equation  given  by  Eq.  5.43, 
Eq.  5.49  may  be  rewritten  as 


U(X, Y) 


3U (X, Y)  3U (X, Y)  fY  3U(X,Y)  . 

jx  5y  oj  Six 


Since  the  jet  is  self-preserving. 


3UV(X,Y) 

3Y 

(5.50) 


.  q(Y/B>  -  ,1,1  (5.51) 

°™<x) 

m 


and  Eq.  5.44  may  be  used  to  write  Eq.  5.50  in  terms  of  the 
mean  velocity  profile  f(n).  After  some  simplification,  the 
result  is 


dq(H)  „ 

d*l 


[BOO 

lum(X) 

m 


dUm(X) 


dX 


dB  (X) 
dX 


]  HP  oi1 


(<})d  q 


B  (X) 

U(X) 

m 


dU  (X) 
m 

dX 


(5.52) 


Integration  of  Eq.  5.52  will  then  yield  the  Reynolds  stress 
distribution  q(n)*  This  procedure  has  been  utilized  by 
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Gutmark  &  Wygnanski  (1976)  with  the  results  given  in 
Figure  5-9.* 

From  Figure  5-9,  it  is  noted  that  the  Reynolds 
stress  calculated  by  the  vortex  street  model  is  in  close 
agreement  with  measurements  and  that  derived  from  the 
Navier-Stokes  equation.  The  capability  of  the  model  to 
duplicate  experimentally  determined  self-preserving  mean 
longitudinal  profiles  and  simultaneously  satisfy  the 
requirements  of  the  continuity  and  Navier-Stokes  equations 
provides  substantial  support  for  the  hypothesized  structure. 
4 .  Velocity  Fluctuation  Intensities 

the  velocity  fluctuation  intensities 

{h  [U (X, Y)  -  U  (X,  Y) ] 2  dt}1/2  (5.53) 


{l  _/*  [V (X, Y)  -  V (X , Y) ] 2  dt}1/2. (5.54) 
T  ^ 

The  resulting  distributions  together  with  experimentally 
determined  ones  are  given  in  Figures  5-10  and  5-11. 

Not  surprisingly,  it  is  in  these  properties  that  the 
vortex  street  model,  as  applied,  fails  to  generate  accurate 
representations  of  the  real  flow.  Although  the  calculated 
distributions  are  qualitatively  reasonable,  the  magnitude 
of  the  lateral  fluctuations  are  about  twice  the  experi¬ 
mental  values  and  the  longitudinal  fluctuations  near  the 
axis  are  much  lower  than  the  experimentally  determined  ones. 


*  The  integration  could  have  also  been  carried  out  with 
the  present  data  giving  similar  results. 


Calculations  of 
were  made  using 

IT(X,Y)  _  1 

0m(x)  "  Sm(X> 

and 

V' (X, Y)  =  1 

U  (X)  U  (X) 
m  m 
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These  problems  are  directly  related  to  the  total 
neglect  of  the  small  scale  turbulence  and  the  simplicity 
of  the  Rankine  vortex  model.  The  extremely  high  values 
of  the  lateral  fluctuation  levels  are  a  result  of  the 
linear  superposition  of  the  velocities  induced  at  the  mid¬ 
point  between  adjacent  counter-rotating  vortices  on 
opposite  sides  of  the  centerline.  In  these  regions,  the 
total  velocities  are  approximately  twice  what  either  vortex 
would  induce  alone.  Undoubtedly,  the  inclusion  of  small 
scale  turbulent  diffusion  effects  would  prohibit  the 
unreasonably  high  local  strain  rates  calculated  by  the 
current  vortex  street  model. 

There  are  two  approaches  that  might  be  utilized  in 
order  to  account  for  the  small  scale  turbulent  effects. 

The  first  is  to  attempt  to  model  the  local  large  scale/ 
small  scale  interaction,  energy  exchange  and  energy 
dissipation  by  the  small  scales,  which  is  not  a  trivial 
task.  The  second  and  more  simplistic  approach  would  be  to 
modify  the  basic  Rankine  vortex  vorticity  and  velocity 
distributions  into  nonaxi symmetric  ones  that  accounted  for 
local  variations  of  the  turbulent  momentum  transport. 

Neither  approach  was  adopted  in  the  currently 
reported  research  since  the  large  scales  were  of  principal 
interest  and  the  Rankine  vortex  was  apparently  sufficient 
to  demonstrate  their  effect  on  mean  velocities  and 
Reynolds  stress  production.  The  key  question  concerning 
the  reported  fluctuation  intensities  is  whether  or  not  they 
are  critical  to  the  apparent  accuracy  of  the  calculated 
U,  V  and  UV  distributions.  This  problem  is  considered  in 
Section  E. 

5 .  Limitations  of  the  Model 

Due  to  the  manner  in  which  the  numerical  simulation 
was  formulated,  several  of  the  experimental  turbulent  flow 
properties  either  cannot  be  calculated  or  their  calculation 
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yields  meaningless  results.  The  primary  missing 
ingredients  in  the  formulation  are  the  effects  of  random 
"jitter"  in  the  coherent  structural  pattern  and  the  effects 
of  the  random  small  scale  turbulent  velocity  fluctuations. 

For  instance,  the  turbulent  kinetic  energy  equation 
cannot  be  evaluated  since  there  is  no  representation  of 
the  small  scales  which  receive  energy  from  the  large  scales 
or  the  even  smaller  scales  which  provide  the  eventual 
viscous  dissipation  of  the  energy.  Similarly,  nothing  may 
be  inferred  concerning  the  turbulent  energy  spectra  since 
the  model  only  represents  one  frequency  component. 

Although  qualitative  agreement  with  the  experimentally 
determined  correlation  map  of  Figure  4-5  is  obvious  for 
the  hypothesized  structure,  the  simulation  does  not  give  a 
realistic  corresponding  pattern.  The  regions  of  positive 
and  negative  correlations  would  have  calculated  correlation 
coefficients  of  +  1  since  the  simulation  is  perfectly 
deterministic. 

While  the  vortex  street  structure  would  determine  the 
overall  interface  geometry,  the  exact  location  cannot  be 
simulated  without  the  small  scale  effects  that  actually 
delineate  the  interface. 

E.  The  Mechanisms  of  Entrainment  and 
Reynolds  Stress  Production 

Considering  that  the  effects  of  the  random  small  scale 
turbulence  were  not  included  in  the  numerical  simulation  of 
the  vortex  street  model,  the  capability  of  the  model  to 
represent  growth  rates,  velocity  decay  rates,  mean  velocity 
profiles  and  the  averaged  Reynolds  stress  profile  is 
remarkable.  It  is  especially  surprising  that  the  calcula¬ 
tion  of  these  properties  was  not  affected  by  the  inability 
of  the  model  to  duplicate  the  experimental  distributions  of 
velocity  fluctuation  intensity.  Browand  &  Weidman  (1976) 
noted  the  invariance  of  mean  flow  properties  and  Reynolds 
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stress  distributions  with  respect  to  Reynolds  number  (and 
the  associated  small  scales)  in  a  plane  mixing  layer.  The 
explanation  which  they  offered  was  that  these  properties 
are  controlled  solely  by  the  large  scales  since  only  they 
are  unaffected  by  the  magnitude  of  the  Reynolds  number. 

The  same  explanation  may  be  offered  for  the  success  of  the 
vortex  street  model  where  it  represents  the  limiting  case 
in  which  the  small  scales  are  totally  absent.  The  question 
then  becomes  one  of  explaining  the  manner  in  which  the 
large  structures  do  generate  the  observed  jet  properties. 

The  relationship  between  the  calculated  mean  flow 
properties  and  the  streamwise  development  of  the  jet  was 
investigated  by  modifying  the  numerical  model.  In  the 
modification,  the  vortex  street  was  represented  as  a 
parallel  flow  with  identical  and  equally  spaced  vortices 
extending  up  and  downstream  to  infinity.  This  is  the 
pattern  that  an  observer  located  at  a  fixed  location  would 
predict  having  only  the  structural  information  available 
at  that  point.  The  calculation  results  duplicated  the 
experimental  mean  longitudinal  velocity  profile  but 
calculated  zero  distributions  of  mean  lateral  velocity  and 
Reynolds  stress.  The  conclusion  drawn  from  this  experi¬ 
ment  is  that  the  global  development  of  the  large 
structures  is  critically  important  to  local  mean  flow 
properties. 

Insight  into  the  effects  of  the  streamwise  development 
of  the  vortex  street  on  entrainment  and  Reynolds  stress 
production  is  gained  from  examination  of  the  instantaneous 
longitudinal  velocity  profiles  given  in  Figure  5-12.  The 
analogous  distributions  for  the  parallel  flow  representa¬ 
tion  are  given  in  Figure  5-13.  The  most  obvious  difference 
between  the  parallel  and  developing  vortex  streets  is  the 
lack  of  symmetry  of  the  distributions  at  t  =  T/4  and 
*t  =  3T/4  for  the  developing  pattern. 
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The  significance  of  the  difference  is  clarified  by 
comparing  the  time  histories  of  the  longitudinal  and 
lateral  velocity  fluctuations  and  Reynolds  stress  at  a 
fixed  point,  (X,Y) .  Plots  of  this  type  are  given  in 
Figure  5-14.  The  fixed  point  (X  =  45,  Y/B  =  0.7)  for  this 
case  corresponds  to  the  lateral  location  where  the  maximum 
Reynolds  stress  production  occurs.  The  key  feature  to  be 
noted  from  the  figure  is  the  slight  decrease  in  the  phase 
lag  between  the  longitudinal  and  lateral  fluctuations. 
Apparently,  the  phase  change  is  a  direct  consequence  of  the 
spreading  rate  and  the  key  to  nonzero  time  averaged 
entrainment  and  Reynolds  stress  production. 

The  dependence  of  the  lateral  velocity  and  Reynolds 
stress  distributions  on  the  spreading  rate  is  implied  by 
the  continuity  and  Navier-Stokes  equations  when  expressed 
in  terms  of  similarity  variables.  From  Eq.  5.46  the 
distribution  of  V(X,Y/B) /Um(X)  is  given  by 
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f(n)  a  U(X,Y/B)/Um(X)  . 


Substituting  Eq.  5.40  for  B(X)  and  Eq.  5.41  for  Um(X) 
yields 

g<n>  -  <J  f<n>  d7  +  (5'56) 

Eq.  5.56  illustrates  that  the  lateral  velocities  are  scaled 
directly  by  the  growth  rate,  B* 


105 


Substituting  for  B(X)  and  Um(X)  in  the  streamwise 
Navier-Stokes  equation  given  by  Eq.  5.52  yields 


q(n)  =  UV(X,Y)/Vm(X)2 


Again,  the  growth  rate  is  noted  to  determine  the  magnitude 
of  the  Reynolds  stress  distribution. 

The  growth  rate,  8,  is  an  experimentally  determined 
coefficient  approximately  equal  to  0.1  for  plane  jet  flows. 
There  has  been  no  purely  theoretical  analysis  reported  which 
was  capable  of  predicting  3.  Since  its  value  was  incor¬ 
porated  into  the  vortex  street  model  as  a  scaling  length, 
it  is  possible  to  use  the  model  as  a  means  of  postulating 
the  flow  fields  which  would  result  for  varying  values  of  8. 
The  only  restriction  is  that  the  accompanying  velocity 
scale  be  compatible.  Conservation  of  momentum  provides  the 
desired  relationship.  Consider 

~  _Jm  U  (X,  Y)  2  dY  =  0  (5.58) 

Since  the  flow  is  self-preserving,  Eq.  5.58  may  be  rewritten 
as 

~  (um(x)2B(x)  _J°°  f (n) 2  dn)  =  0 
or 

U  2 (X) B (X)  =  constant  (5.59) 

m 
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Substituting  Eqs.  5.41  and  5.42  for  B  and  Um»  respectively, 
yields 

C^B  =  (2.36)2  *  0.1  (5.60) 

where  the  constant  has  been  determined  from  the  naturally 
occurring  flow. 

Figures  5-15  through  5-17  illustrate  distributions  of 
U,  V  and  UV  calculated  by  the  vortex  street  model  for  a 
wide  range  of  jet  growth  rates.  For  comparison,  the 
corresponding  distributions  calculated  from  Eqs.  5.56  and 
5.57  are  also  given. 

From  data  similar  to  that  shown  in  Figure  5-14,  the 
phase  change,  6$,  was  determined  for  each  widening  rate 
with  the  results  given  in  Figure  5-18.  Surprisingly,  6<f> 
and  6  are  approximately  linear  with  respect  to  one 
another. 

In  a  remarkably  insightful  publication,  Townsend 
(1979)  has  postulated  that  the  global  flow  development 
introduces  the  same  type  of  phase  shift  effect  between  the 
velocity  components  in  a  two-dimensional  wake.  He  noted 
that  the  velocity  product  of  fluctuations  represented  by 
sinusoids  exactly  180°  out  of  phase  would  not  on  the 
average  generate  a  nonzero  Reynolds  stress.  If  from  the 
overall  flow  development,  a  slight  phase  change  was 
introduced,  a  nonzero  product  would  result.  The  present 
results  substantiate  that  effect. 

The  distributions  of  the  fluctuating  velocity 
intensities  for  the  various  widening  rates  are  given  in 
Figures  5-19  and  5-20.  It  is  noted  that  they  are  scaled 
by  Tjm  and  are  not  dependent  on  B.  Consequently,  it  is 
apparent  that  the  capability  of  the  vortex  street  model  to 
generate  realistic  distributions  of  entrainment  velocity 
and  Reynolds  stress  is  not  dependent  on  the  corresponding 
fluctuation  intensity  distributions. 
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The  fluctuation  intensity,  mean  velocity,  and  Reynolds 
stress  profiles  are  related,  however.  For  a  given  approxi¬ 
mate  model  of  the  individual  vortices  making  up  the  vortex 
street,  there  is  a  corresponding  set  of  similarity  con¬ 
stants,  CH,  CR  and  Cp,  which  optimize  the  agreement  with  the 
mean  longitudinal  velocity  profile  (C_  is  fixed  from  experi- 
mental  observations) .  The  accuracy  with  which  the  corres¬ 
ponding  Reynolds  stress  and  lateral  velocity  profiles  may 
be  calculated  is  apparently  directly  related  to  the 
accuracy  with  which  the  associated  mean  velocity  profile  is 
duplicated  and  independent  of  the  associated  fluctuation 
intensities.  The  inaccuracy  of  the  fluctuation  intensities 
is  therefore  an  unavoidable  side  effect  of  utilizing  a 
simplistic  vortex  model  such  as  the  axisymmetric  Rankine 
vortex. 

G.  Summary 

A  numerical  simulation  of  the  vortex  street  model  for 
the  large  scale  structures  of  a  turbulent  plane  jet  has  been 
described.  The  formulation  includes  the  effects  of  stream- 
wise  flow  development  through  the  expression  of  geometry 
and  vortex  parameters  in  terms  of  similarity  length  and 
velocity  scales. 

The  calculation  results  indicate  that  the  hypothesized 
structure  is  capable  of  giving  virtually  exact  representa¬ 
tions  of  the  growth  rate,  velocity  decay  rate,  mean  velocity 
profiles  and  Reynolds  stress  distributions  of  real  jet 
flows.  Although  the  qualitative  shapes  of  the  fluctuation 
intensity  profiles  were  in  agreement  with  experimentally 
determined  ones,  the  magnitudes  were  considerably  in  error. 
It  was  argued  that  this  was  a  consequence  of  the  failure  to 
include  the  effects  of  the  small  scale  turbulence. 

A  key  feature  of  the  combined  calculated  mean  velocity 
and  Reynolds  stress  profiles  was  that  taken  together,  they 
demonstrated  compatibility  with  the  requirements  of  the 
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two-dimensional  continuity  and  the  Navier-Stokes 
equations.  This  capability  is  particularly  notable 
considering  the  simple  kinematic  representation  of  the 
flow  field. 

Detailed  examination  of  the  relationship  between 
instantaneous  Reynolds  stress  and  the  velocity  fluctuations 
revealed  that  the  downstream  development  of  the  jet  flow  is 
critical  to  the  local  entrainment  and  Reynolds  stress 
production  rates.  This  dependence  on  the  growth  rate  is 
implicit  in  the  continuity  and  Navier-Stokes  equations. 

The  model  successfully  represents  cases  in  which  the  growth 
rate  is  different  from  that  occurring  naturally  in  plane 
jets.  The  strongest  outward  indication  of  the  effect  of  a 
nonzero  growth  rate  is  the  introduction  of  a  change  in  the 
phase  relationship  between  the  longitudinal  and  lateral 
velocity  fluctuations. 

In  Section  V.A  ,  it  was  noted  that  an  alternative 
representation  of  the  vortex  model  might  have  been  the 
superposition  of  the  structure  onto  an  otherwise  homo¬ 
genous  mean  flow.  From  the  preceding  analysis,  it  is 
apparent  that  this  would  have  upset  the  balance  between  the 
longitudinal  mean  velocity  profile  and  the  requirements  of 
continuity  and  the  Navier-Stokes  equations.  It  is 
necessary  for  the  large  scale  structures  to  possess 
sufficient  energy  to  duplicate  the  mean  longitudinal 
velocity  profile  before  the  appropriate  accompanying 
lateral  velocity  and  Reynolds  stress  distributions  can  be 
generated. 
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CHAPTER  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  primary  focus  of  the  research  which  has  been 
reported  was  to  gain  insight  into  the  characteristics  and 
importance  of  the  large  scale  structures  in  the  similarity 
region  of  a  two-dimensional  turbulent  jet.  This  was 
motivated  by  experimental  indications  early  in  the  program 
of  a  vortex  street  type  of  coherent  eddying  pattern.  The 
vortex  street  model  was  then  adopted  as  a  working 
hypothesis  for  the  interpretation  of  the  experimental  and 
numerical  analysis  that  followed. 

A.  Conclusions  Based  on  Experimental  Results 
Several  conclusions  relative  to  the  hypothesized 
structure  may  be  drawn  from  the  experimental  results  which 
have  been  reported. 

1.  Every  experimental  indication  points  to  the  fact  that 
the  concept  of  structural  similarity  is  applicable  to 
the  large  structures  of  a  fully  developed  jet. 

Support  for  this  conclusion  is  drawn  from  the 
similarity  scaling  of  the  intermittent  region  and  the 
integral  length  scales.  In  addition,  the  streamwise 
variation  of  the  interface  crossing  frequency  is 
compatible  with  similarity  considerations. 

2.  Comparison  of  the  interface  crossing  frequency 
(measured  using  a  rationally  defined  hold  time 
variation)  with  the  previously  reported  "flapping" 
frequencies  indicates  that  the  interface  geometry  and 
underlying  structure  are  related.  Further  comparisons 
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with  reported  "flapping"  frequency  measurements 
indicate  a  structural  frequency  constant, 
fsb/Um  =  0.1  for  all  plane  jet  flows.  This  suggests  a 
universal  instability  mode  which  leads  to  the  vortex 
street  structure. 

3.  Two  point  correlations  within  the  turbulent  velocity 
field  yield  the  same  spatial  relationships  which  would 
be  expected  of  a  vortex  street-like  structure.  These 
include  a  characteristically  antisymmetric  instanta¬ 
neous  longitudinal  velocity  profile,  streamwise 
periodicity  and  coherence  over  substantial  distances. 

It  is  not  clear  how  these  patterns  develop  within  the 
initial  and  transition  regions  (i.e.,  0  <  x/D  <  10). 
This  is  one  area  of  suggested  further  research. 

4.  Utilizing  similarity  concepts  and  relationships  based 
upon  the  premise  of  a  vortex-like  structure,  the 
convective  velocity  of  the  structural  pattern  was 
determined  from  flow  visualization  experiments.  The 
result  is  that  the  pattern  (i.e.,  the  vortex  centers) 
travels  at  approximately  one-half  the  local  mean 
centerline  velocity. 

5.  Efforts  at  determining  the  extent  of  two-dimensionality 
characteristic  of  the  large  vortex-like  structures  led 
to  inconclusive  results.  Comparison  of  axial  and 
vertical  integral  length  scales  within  the  intermittent 
region  suggest  at  least  limited  two-dimensionality. 

The  same  is  not  true  for  the  integral  scales  of  the 
turbulent  velocity  field.  There  is,  however,  an 
indication  from  the  vertical  correlation  measurements 
that  a  vortex-like  structure  may  exist  which  is  more 
complex  than  the  hypothesized  structure,  yet  still 
substantially  coherent  in  the  vertical  direction.  It 
is  suggested  that  this  feature  be  investigated  further 
possibly  utilizing  numerically  applied  pattern 
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recognition  techniques  such  as  those  used  by 
Blackwelder  (1977)  in  the  two-  dimensional 
turbulent  boundary  layer. 

B.  Conclusions  Based  on  the  Numerical  Simulation 

Based  upon  the  results  of  the  experimental  investiga¬ 
tion,  a  numerical  simulation  of  a  two-dimensional  jet  formed 
from  a  vortex  street  was  evaluated.  The  following 
conclusions  were  reached: 

1.  A  vortex  street  model  is  quite  capable  of  duplicating 
experimentally  determined  mean  velocity  profiles, 
growth  rates,  velocity  decay  rates  and  Reynolds  stress 
distributions  but  fails  to  calculate  intensities, 
spectra  and  correlations. 

2.  Flow  properties  calculated  by  the  numerical  simulation 
satisfy  the  two-dimensional  approximations  of  the 
fundamental  conservation  requirements  expressed  by  the 
continuity  and  Navier-Stokes  equations. 

3.  It  is  evident  that  the  key  feature  of  the  vortex 
street  model,  as  formulated,  is  the  inclusion  of  the 
effects  of  global  flow  development  on  locally 
measured  or  calculated  flow  properties. 

4.  As  a  consequence  of  the  flow  development,  a  change  in 
the  phase  relationship  between  longitudinal  and 
lateral  velocity  components  is  introduced.  It  is 
suggested  that  this  feature  be  verified  experimentally 
in  plane  jet  flows  which  have  different  spreading 
rates  due  to  acoustic  forcing  or  some  other  form  of 
external  excitation. 

5.  The  numerical  simulation  suggests  that  increasingly 
negative  longitudinal  velocity  components  at  the 

jet  edges  are  a  consequence  of  the  finite  length 

of  streamwise  structural  coherence.  This  could 
possibly  be  verified  by  examining  the  effect  on  the 
outer  flow  of  artificially  disrupting  the  structure  at 
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increasingly  shorter  downstream  lengths. 

6.  The  effects  of  smaller  random  scales  of  turbulence 
were  not  included  in  the  simulation  and  unrealistic 
distributions  of  flow  properties  directly  related  to 
them  are  the  consequence.  It  is  suggested  that  the 
vortex  street  simulation  serve  as  a  starting  point 
for  the  modeling  of  the  large  scale/small  scale 
interactions  and  turbulent-nontr''1  ulent  interface. 

7.  The  most  astonishing  implication  of  the  model  is  that 
entrainment  and  Reynolds  stres?  production,  which  are 
key  elements  of  any  free  turbulent  shear  flow,  are 

i 

actually  consequences  of  the  coherent  or  nonturbulent 
aspects  of  the  flow.  This  is  actually  a  question  of 
semantics  which  serves  to  illustrate  the  type  of 
redefinition  of  fundamental  turbulence  concepts  which 
must  result  from  the  increasingly  convincing  indica¬ 
tions  of  the  importance  of  coherent  structures  in 
turbulence. 

8.  As  noted  on  page  95,  the  model  predicts  a  negative 
distribution  of  Reynolds  stress  although  a  positive 
distribution  is  required  by  the  Navier-Stokes  equation. 
This  problem  is  clearly  indicated  in  Figure  5-14.  A 
thorough  search  has  been  made  for  a  simple  sign  error 
in  the  algorithm  but  none  was  found.  The  possibility 
then  arises  that  the  hypothesis  of  a  vortex  street 
structure  in  plane  jets  is  invalid.  This  is  difficult 
to  accept  due  to  the  experimental  evidence  and  the 
capability  of  the  model  to  calculate  U,  V  and  |uv|. 

It  is  suggested  that  there  are  some  subtle  aspects  of 
the  mean  flow  effects  and  energy  transport  processes 
which  have  not  been  properly  accounted  for  in  the 
model  and  that  further  research  be  made  in  this  area. 
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APPENDIX  A 

LITERATURE  REVIEW  -  TURBULENCE  DETECTION 


In  principle,  the  formulation  of  the  intermittency 
function  and  the  measurement  of  the  intermittency  fraction 
and  interface  crossing  frequency  are  straightforward.  How¬ 
ever,  in  practice,  the  procedure  is  difficult  due  to  the 
sensitivity  of  I(r,t)  to  the  choices  made  for  the  original 
indicator  signal,  threshold  and  hold  time.  Most  variations 
on  Townsend's  original  intermittency  circuit  have  been 
conceived  with  the  intent  of  overcoming  these  difficulties. 
What  follows  is  a  summary  of  the  techniques  that  have  been 
utilized  and  reported  in  the  literature. 


A.  The  Intermittency  Indicator  Signal,  S(r,t) 
Different  choices  for  S(r,t)  have  been  tried  by 
several  investigators  for  the  purpose  of  clarifying  the 
determination  of  and  C  by  increasing  the  turbulent-non- 
turbulent  contrast  of  S(r,t).  A  representative  sample  of 
these  different  types  of  signals  are  tabulated  later  in  the 
chapter  in  Table  1. 

Townsend's  choice  for  S(r,t)  was  the  rectified  time 


derivative  of  the  streamwise  component  of  the  velocity 
fluctuation.  As  would  be  expected,  the  fluctuation 
intensities  are  appreciably  greater  in  turbulent  than  in 
nonturbulent  regions.  The  contrast  is  enhanced  by  taking 
the  time  derivative  which  favors  the  high  frequency 
fluctuation  components  characteristic  of  the  turbulent 


regions.  This  type  of  signal  has  been  utilized  by  Corrsin 
and  Kistler  (1955),  Fiedler  and  Head  (1966),  Wygnanski  and 
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Fiedler  (1969)  ,  Thomas  (1973) ,  Antonia  and  Atkinson  (1974) 
and  Ajagu  and  Goldschmidt  (1976) . 

The  rectified  time  derivative  of  the  instantaneous  uv 
Reynolds'  stress  component  was  found  to  give  a  more 
distinct  turbulent-nonturbulent  contrast  by  Antonia  (1972) . 

As  suggested  by  Corrsin  and  Kistler,  the  seemingly 
ideal  choice  for  S(r,t)  is  the  instantaneous  vorticity. 

With  the  exception  of  recent  advances  made  by  Foss  (1977) , 
the  practical  difficulties  of  constructing  a  probe  to  sense 
even  a  single  component  of  the  vorticity  vector  have  been 
prohibitive.  Kibens  and  Kovasznay  (1970)  have,  however, 
utilized  a  signal  that  is  approximately  proportional  to  the 
transverse  vorticity  component  by  forming  6u/6y  from  the 
differential  output  of  two  parallel  probes  separated  by  a 
small  transverse  distance.  This  signal  was  differentiated 
with  respect  to  time  and  rectified  to  form  S(r,t).  Jenkins 
and  Goldschmidt  (1974)  have  also  used  this  procedure. 

There  is  an  inherent  difficulty  in  the  use  of  the 
velocity  fluctuation,  its  spatial  or  time  derivatives  or  any 
similar  flow  property  which  fluctuates  between  positive  and 
negative  values  as  a  basis  for  S(r,t).  This  difficulty  is 
due  to  the  frequent  threshold  crossings  that  are  not 
related  to  an  interface  passage  but  are  instead  a  result  of 
the  fluctuating  signal  crossing  zero.  The  occurrence  of 
these  extraneous  threshold  crossings  increases  the  difficulty 
in  choosing  xh  and  adds  to  the  uncertainty  of  intermittency 
measurement. 

To  avoid  the  zero  crossing  problem,  Kaplan  and  Laufer 
(1968)  defined  S(r,t)  as  the  variance  of  the  time  derivative 
of  the  velocity  fluctuation  from  its  short  time  mean  value. 
The  time  interval  over  which  the  mean  value  was  formed  was 
chosen  to  be  small  compared  to  the  expected  burst  and  void 
lengths  and  centered  over  the  instant  of  interest  so  that 
the  past  and  future  were  weighted  equally. 
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Wygnanski  and  Fiedler  (1970)  found  that  the  number  of 
non-interface  related  threshold  crossings  could  be  reduced 
by  forming  S(r,t)  from  the  sums  of  the  rectified  first  and 
second  derivatives  of  the  velocity  fluctuation. 

A  similar  effect  was  achieved  by  Fiedler  and  Head 
(1966)  who  used  the  sum  of  the  rectified  time  derivative  of 
the  velocity  fluctuation  signal  phase  shifted  by  a  time 
interval  that  was  small  compared  to  burst  and  void  durations 
with  the  original  unshifted  but  rectified  signal.  Hedley 
and  Keffer  (1974)  approached  this  problem  by  forming  S(r,t) 
from  the  sum  of  the  squares  of  time  derivatives  of  the 
longitudinal  and  transverse  components  of  the  fluctuating 
velocity  vector. 

A  different  approach  to  the  "zero  crossing"  problem  is 
to  base  S(r,t)  on  the  detection  of  a  marked  species  trans¬ 
ported  within  the  turbulent  fluid  so  that  the  indicator 
signal  will  contain  only  positive  deviations  from  the 
reference  level.  The  advantage  of  having  only  positive 
deviations  is  obvious.  The  disadvantage  is  the  uncertainty 
as  to  the  contamination  of  the  flow  field  by  the  tracked 
scalar.  Freymuth  and  Uberoi  (1973) ,  LaRue  (1974) ,  Barsoum, 
Hawaii  and  Keffer  (1978)  and  Jenkins  and  Goldschmidt  (1974) 
have  utilized  the  measurement  of  temperature  in  slightly 
heated  turbulent  flows  for  this  purpose.  A  helium  marked  jet 
has  been  investigated  by  Anderson,  LaRue  and  Libby  (1977) 
while  Fiedler  and  Head  (1966)  revealed  the  intermittent 
characteristics  of  a  smoke  filled  boundary  layer  by 
illuminating  it  and  observing  the  alternating  light  and  dark 
patches  with  a  high  resolution  optical  system. 


B.  The  Intermittency  Detection  Parameters,  C  and 

The  greatest  uncertainty  in  forming  I(r,t)  is  in 
defining  the  criteria  for  selecting  the  threshold*,  C,  and 
hold  time  t^.  Most  investigators  have  achieved  satisfactory 
results  by  considering  each  parameter  separately. 

The  optimum  choice  for  the  threshold  is  the  minimum 
value  which  is  stilll above  the  inherent  "noise"  contained 
in  the  nonturbulent  portions  of  S(r,t).  Wygnanski  and 
Fiedler  (1969),  Antonia  (1972),  Kaplan  and  Laufer  (1968), 
and  LaRue  (1974) ,  have  found  that  this  level  was  most 
easily  determined  through  subjective  inspections  of  S(r,t) 
and  the  I(r,t)  which  resulted  from  a  particular  choice  of  C 
as  shown  by  oscilloscope  or  computer  plotted  traces. 

Jenkins  and  Goldschmidt  (1974)  ,  and  Ajagu  and  Goldschmidt 
(1976) ,  determined  C  by  observing  the  variation  y  and 
for  a  range  of  thresholds.  Although  y  systematically 
decreased  for  increasing  C,  the  rate  of  variation  of  F^ 
decreased  markedly  for  thresholds  above  a  value  which  was 
interpreted  as  the  optimum  value.  A  similar  approach  was 
taken  by  Hedler  and  Keffer  (1974)  using  the  variation  of  the 
average  nonturbulent  time  interval  as  the  critical 
threshold  dependent  property. 

Kibens  and  Kovasznay  (1970) ,  found  that  the  optimum 
choice  for  the  intermittency  of  I(r,t)  was  a  function  of  the 
intermittency  and  the  ratio  of  the  zone  averaged  rms  levels 
of  the  turbulent  and  nonturbulent  portions  of  S(r,t)  which 
was  defined  as  the  signal  to  noise  ratio.  Since  the  inter¬ 
mittency  and  signal  to  noise  ratio  are  themselves  dependent 
on  C,  calibration  plots  for  the  optimum  C  were  determined 

*  The  threshold  is  typically  defined  as  either  the  percentage 
of  some  characteristic  measured  flow  property  (peak  or  rms) 
or  alternately  as  an  absolute  voltage  level.  The  hold  time 
is  always  given  in  units  of  time. 
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from  analysis  of  an  artificial  turbulence  signal  with  known 
intermittency  properties.  An  iteration  procedure  could  be 
used  to  find  C  by  assuming  a  value  of  y  and  the  S/N  ratio 
and  picking  the  optimum  threshold  from  the  calibration 
plots.  The  measured  y  and  S/N  were  compared  to  the  assumed 
values  and  additional  iterations  made  accordingly. 

The  hold  time  is  typically  determined  from  estimations 
of  the  minimum  expected  length  scale,  ^  ,  of  the  eddy 
structure  associated  with  interface  deformations,  is 

then  approximately  equal  to  ^min/uc  where  Uc  is  the  convec¬ 
tion  velocity.  Hedley  and  Keffer  (1974) ,  and  LaRue  (1974) , 
expressed  as  an  integer  multiple  of  the  Kolmogoroff 

length,  Kibens  and  Kovasznay  (1970) ,  Jenkins  and 

Goldschmidt  (1974) ,  and  Antonia  and  Atkinson  (1954) ,  found 
that  could  be  expressed  in  terms  of  a  characteristic 
frequency  defined  as  F  ,  =  S (r , t) /S (r , t) . 

A  completely  different  approach  was  adopted  by  Ajagu 

and  Goldschmidt  (1976) ,  who  chose  t.  as  the  minimum  value 

h 

that  would  allow  measurement  of  an  intermittency  equal  to 
one  in  a  fully  turbulent  zone  for  the  range  of  thresholds 
which  were  used.  Unlike  other  investigators  who  based 
on  the  measurement  of  local  properties  so  that  it  had  a 
nonconstant  value,  Ajagu  and  Goldschmidt  (1976)  were  able 
to  use  a  constant  over  a  wide  portion  of  the  self¬ 
preserving  region  of  a  plane  jet. 

The  most  involved  procedure  for  the  determination  of 

and  C  was  developed  by  Antonia  and  Atkinson  (1974) .  They 

considered  the  optimum  choice  of  and  C  to  be  functionally 

related  to  each  other  as  well  as  to  the  intermittency  and 

S/N  ratio.  Using  calibration  plots  derived  from  the 

analysis  of  artificial  turbulent  signals,  C  and  t.  were 

h 

determined  through  an  iteration  procedure  similar  to  that  of 
Kibens  and  Kovasznay. 


A  partial  tabulation  of  a  few  of  the  intermittency 
measurement  procedures  is  given  in  Table  A-l.  The 
uncertainty  of  the  measurement  process  is  attested  to  by 
the  wide  variety  of  techniques. 


\ 


determined  from  subjective  comparisons  of  S(r,t)  and  resulting  I(r,t) 
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APPENDIX  B 

THE  DIGITAL  INTERMITTENCY  MEASUREMENT  TECHNIQUE 

Although  the  plane  jet  apparatus  is  located  in  the 
R.  W.  Herrick  Laboratories,  the  analysis  of  the  early 
intermittency  data  was  done  on  a  General  Automation  mini¬ 
computer  located  in  the  mechanical  engineering  stochastics 
laboratory.  For  this  reason,  the  raw  turbulence  indicator 
signal  was  recorded  on  an  Ampex  FM  tape  recorder  and  later 
played  back  for  digitization  and  analysis.  A  block  diagram 
of  the  data  acquisition  and  analysis  equipment  is  given  in 
Figure  B-l. 

Analysis  of  the  digitally  recorded  data  was  accom¬ 
plished  by  subroutine  DTECTR  under  the  control  of  the  main 
routine  TDAP .  Listings  of  these  routines  are  given  at  the 
end  of  this  section.  The  logic  for  the  application  of  the 
threshold  (designated  as  LEVEL  in  DTECTR)  and  hold  time 
(designated  as  HOLDTM)  may  be  given  in  three  rules: 

1.  Data  samples  ( DATA ( 1 ) )  greater  than  LEVEL  are 
identified  as  being  part  of  turbulent  bursts  while 
samples  below  LEVEL  represent  voids  or  the  intervals 
separating  bursts. 

2.  All  bursts  and  voids  must  have  durations  greater  than 
the  hold  time  to  be  valid,  i.e.,  a  burst  must  last 
longer  than  the  hold  time  before  it  may  interrupt  an 
otherwise  continuous  void  and  vice  versa. 

3.  A  conflicting  case  may  arise  for  rule  #2  when  a  series 
of  short  bursts  are  separated  by  short  voids  with 
neither  the  bursts  nor  voids  longer  than  the  hold  time. 
For  this  case,  the  short  bursts  and  voids  are  counted 
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as  a  single  burst  provided  the  total  duration  is  then 
greater  than  the  hold  time. 

A  flow  chart  describing  the  implementation  of  these  rules 
and  an  example  of  their  application  is  given  in  Figure  3-3. 


i 
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Table  B-l  Program  TDAP 


TDAP  -TURB.  DATA  ANALYSIS  PROGRAM  1  4  M  F  J*0 

DIMENSION  y  C  3wJ )  .  FREQ130).  HE  I  GUT  (  2 )  .  T  HRHLD l  30) 

REAL  INTERM(30).  MXFREu  .  NSET 

INTEGER  SIRTVE.  NFILE130).  OPTION 

INTEGER  HEADER.  LABt  A  ( 4  ) .  LAB.  *1(8).  LAB.  Y2 l 1 0) 

DATA  HEADER  70/ 

DATA  LAfcLX  /  *Y  * . !/  * . ‘B  ‘ .2/ 

DATA  LAbLYl  /  *  I N ‘ . * TE * . * RM * .  *  I T * . ' TE * . * NC 1 .  *  1  ’.1/ 

DATA  LAbi  il  /  *CR ' . 'OS ’ . *  SI  * . *NG ‘ *  F * . • Rt * . * QU  ’  .  ’ EN 1 . ' CY 

09  CONTINUE 

ARITEl 7. 107) 

READ! 7. 999)  NFILES 
IFINriLES.EQ.0)GO  TO  98 
ARI TE(7> 101 ) 

READ ( 7.999)  HOLDTM 
JYRlTEl  7.  Ill) 

READ (7.999)  DT 
#RHEl7.  108) 

READ! 7.999)  NSET 
NN=NFILES+2 
A  R I T  E  (  7 .  109) 

DO  9  1=3. NN 

J=  I  -2 

#RITE<  7. 110)  J 

READ (  7.9  99)  Y (I)  ,NFILE(I) , THRHLD ( I ) 

9  CONTINUE 

CAL-ULATE  IN  I ERM  I TTENCY  fi,  CROSSING  FREQ.  DISTRIBUTIONS 
DO  10  1=3. NN 

CALL  DiOPEN 

ILENGH  =  STRTVEtNFlLE(I). 14.0. 0.-1) 

CALL  DiOPEN 

CALL  DTEC I  RIOT .HOLDTM. THRHLD ( I ) . NSET . INTERM ( I  )  . FREQ ( I ) ) 

10  CONTINUE 

NORMALIZE  ThE  CROSSING  FREQ.  DISTRIBUTION 
MXFREQ  3  0.0 
DO  11  1=3. NN 

IFlFREu(I)  .GT.  MXFREQ)  MXFREQ  =•  FREQ ( I ) 

11  CONTINUE 

DO  12  I  =  3.NN 

FREQ!  1  )  ••  FREQ  1  l  )  /  MXFREQ 

12  CONTINUE 
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AfRITE  IHE  RtSULlS 
H/R1TE15. 102) 

lRIfEt5< 103)  01.  HOLDTM.  MXf-REu 

#R1TE!5. 104) 

*R1TE15.105)  !  Y!  I)  .  THRHLDU  >  .  INTERMU).  FREQ  Cl).  1*3.  NT 

PLOT  THE  Rl SJLTS 
«R1TL(7.  10-3) 

READ (7.990)  OPTION 

IF  <  OPTION  .LT.  0  )  GO  TO  99 

IF  l  OPTION  ,NE.  1  >  GO  TO  13 

Y! 1 )=>0.0 

Y!2)=2.0 

INTERVll)  *  0.0 

1 NTERM (2 )  *  1.0 

FREQ( 1)  *  0.0 

F REQ ( 2 )  *  1.0 

LAbLXil)  *  *  Y  ' 

LABLYl(l)  *  ‘IN* 

HEIGHT ( 1 )  *  5.0 
HEIGHT ( 2 )  =  5.0 

CALL  JPLOT 1  1  .NFILES. -1 . 1 . HEADER . Y . LABLX , -959 . 0 . I NTERM .L A£ 
X  -999. 0, HEIGHT .-1 .0) 

LABLY2U)  *  *  CR  ‘ 

HEIGHT! 1 )  *  5.0 
HEIGHT  C2 )  s  0.0 

CALL  JPLOT  10.  NFILES.  -  1  .  1 .  HEADER  .  Y  .  LABLX  .  -99‘J  » 0  .  FREQ  .  LABL\ 
*  -999. 0. HEIGHT. -1  .0) 

GO  TO  99 

13  CONT  I.nUE 

HEIGHT! 1 )  =  5.0 
HEIGHT ! 2 )  *  5.0 
LAtJLXIl)  0 
LAbLYl ! 1 )  =  0 

CALL  JPLOT  t  0.NEILES.-1  .  1 .  HEADER  .  Y  .LABL.X .  999 . 0 .  INTERM .  LAbL 
A  HEIGHT .-1.0) 

HEIGHT! 1 )  5.0 

HEIGHT  12 )  ^  0.0 
LAULY2! 1 )  *  0 

CALL  JPLOT l 0. NFILES. -1 . 1 .HEALER . Y. LABLX . 990. 0. FREQ. L ABLY 
X  HEIGHT. -1.0) 


GO  To  99 


95  CONTINUE 


4  V 
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Table  B-l  cont. 


101  format!'  HOLUTM?*) 

102  FORMAT!///'  INJERMITTENCY  &  CROSSING  FREQ  •  DISTRIBUTIONS' 
10J  FORMAT !  /4X .  *01  ^  '  .El  1 *3.4X. ‘HOLD  TIME  -  *.F11.3./4X. 

X ' MAX  CROSSING  FREQ  =  *.E11.3//> 

104  F0RMATI//9X. * Y/B * , 8X, *  THRHLO ' ,5X, '  INTERM  I TTENCY * ,3X, *CR. 

105  F0RMAT14 ( JX.El 1 .3) > 


106  FORMAT!//' 

1 

-  PLOT  ON  A  NEW 

PAGE  '/. 

X  '  ENTER ! 

0 

-  PLOT  ON  SAME 

PAGE'  /. 

X  ' 

107  FORMAT! '  Nf ILLS? ' ) 
106  FORMAT! '  NSET?  ' ) 

-1 

-  CONTINUE  w/o 

PLOT  '  ) 

109  FORMAT (•  Y (I)  ,  NFILE(I), 

THRHLD (I)  ') 

110  FORMAT! IX. 13) 

111  FORMAT!'  OT?  ') 
5)99  FORMAT!  V) 

END 
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Table  B-2  Program  DTECTR 


**  DltCirt  -TURD  DJRSI  DETEC I  ION  ROUT I  Nt  1  4  S  F  J*0  R  PRlE 
SUgROUT 1NE  U I ED 'RIOT .  HOLD! M . LEVEL . NStT . I NTERM . FRLG ) 

0  1  Me RSI  ON  DmTA(J) 

REAL  LEVEL  .JNTtRto  .  NSLT .  ILL.  ITE.  I 
i  =  l 
I  *  0.0 
F  URB  T  m  -  0.  0 


FIND  THE  FIRST  T URQ  bJRS f 

1  CONTINUE 

l  -  I  ♦  1.0 

IFI  1  .bf.  N3ET  )  GO  TO  9131 
CALL  VSREAUDaTA. 1ER) 

IFliER  . NL .  13)  GO  To  950 
1F(DaTA(1)  .GT.  LEVEL)  GO  TO  1 
I LE  =  1 

FIND  THE  TRAILING  EDGE 

2  CONTINUE 

1  3  I  ♦  1 *0 

IK  I  .GT.  NSET )  GO  TO  902 
CALL  VSREALlDATA. ItR) 

1FUER  .NL.  0)  GO  TO  950 
1 F  <  0 ATA (  1  )  .LE.  LEVEL)  GO  TO  2 
1  TE  -  l 

LOOK  FOR  I  HE  Ni  AT  BURST  AND  VER I F f  THAT  INTERVENING  VOID 
IS  LONGER  THAN  HOLDTm 

3  CONI INJE 

i  *  1  ♦  1  •  0 

IF (I  .GT.  NSET)  GO  TO  90S 
CALL  ViREALlDATA. ItR) 

IFIILR  .Nt.  13)  uO  Tu  950 
IF  (DATA!  1  )  .GT.  LEVEL)  GO  TO  3 
VOID  Tm  =  (  1  -  ITL)  *  DT 
IF  CVOIDTM  .LF.  HOlDTM)  GO  TO  2 

VERIFY  THAT  PREVIOUS  T URB  bJRSF  #AS  LONGER  THAN  HOLDTM 
BuRbTM  -  (IVE  -  ILL)  *  DT 
IMiiuRSTM  .LT.  HOLDTM)  GO  TO  4 

1 NCRLMEnI  J  -  INDEX  AND  At D  BURST  LENGTH  FO  (OTAL 
J  -  J  *  1 

lURUTM  -  I  URU I  V  UuRSTm 

BU.lN  CAL  u:  A I  ION  FOR  NE  4  bJRST  PROSPECT 
t  CONI  i,  JE 
I  LE  i 
G 0  ft;  2 
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Table  B-2  cont. 

COMPLETE  Si  ARCH 
00  i  COM  I.'.Jt 
J  •  J 
GO  f  i  j  1  0 

vj «■) a  continue 

1  TE  -  I 
J33  CONTINUE 

TURDTM  -  lUKOIM  *•  CITE  -  ILE)  *  DT 
J10  CONTI’. UE 

C ALviUL  ATE  1NTEHM  I  TANCT  AND  CROSSING  FREQ 
lNTERM  -  T U R 0 T r j*.  /(NSET  *  DT  ) 

FHEO  »  J  /  (NSET  *  DT) 

«RI 1 E  ThE  RlSULTS 

#RI fEl 7. IU1 )  LEVEL. INTERM. FREQ. J 
Ul  FORtaATC  THRHLD  *■  '.  F6 . 3 . 3X  .'  INTtRM  *  *  .  E  1  0. 3 . 3X .  *CR  FRE 
EIJ.3. 3C.  1  J  -  *  .  IS) 

GO  TO  060 

1 02  FORMATC  DAS  INPUT  ERROR  -  I ER  -  *.U) 

9 0  0  CON  I  I  RUE 

#RI I E ( 5 .  102)  IER 
OoU  CONI  I NU £ 

RETURN 
E  NO 
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APPENDIX  C 

THE  HYBRID  INTERMITTENCY  INSTRUMENT 

In  an  effort  to  overcome  the  inconvenience  of 
performing  data  acquisition  and  data  analysis  at  two 

different  locations,  an  alternate  means  of  formulating  the 

•  •  • 

mtermittency  function,  I(r,t),  was  devised.  The  resulting 
hybrid  instrument  has  the  virtues  of  low  cost  and 
convenience  inherent  with  analog  instrumentation  while 
retaining  much  of  the  precision  of  a  larger  digital  computer. 
To  the  author's  best  knowledge,  the  instrument  represents 
the  first  application  of  a  computer  of  any  size  for  the 
continuous  and  real-time  formulation  of  the  intermittency 
function. 

1 .  The  Hardware 

A  block  diagram  of  the  basic  elements  of  the  inter¬ 
mittency  detection  instrument  is  given  in  Figure  C-l.  They 
consist  of  a  hot  wire  anemometer,  DC  signal  conditioning 
circuits,  the  threshold  circuit  and  an  RCA  Cosmac  micro¬ 
computer. 

The  DC  signal  conditioning  circuits  consist  of  a 
reverse  biasing  circuit  to  remove  the  DC  component  of 
S(r,t)  and  a  DC  amplifier.  Application  of  C  to  S(r,t)  is 
accomplished  with  an  operational  amplifier  comparator.  With 
S(r,t)  at  the  negative  input  and  a  variable  C  at  the  positive 
input,  a  square  wave,  G{r,t)  is  available  at  the  output. 
S(r,t),  C  and  G{r,t)  are  all  available  as  outputs  to  an 
oscilloscope  which  facilitates  the  adjustment  of  C. 


■7*- 
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Figure  C-l  The  Hybrid  Intermittency  Analysis  Technique 
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The  microcomputer  is  based  on  the  RCA  Cosmac  family 
of  microprocessor  integrated  circuits  using  an  eight  bit 
data  bus  and  16  bit  address  bus.  Operating  system  firmware 
is  contained  in  1/4K  ROM  and  3K  RAM  is  available  for  user 
programs.  An  audio  cassette  interface  is  utilized  for 
program  storage  and  retrieval.  The  intermittency  circuit  is 
interfaced  to  the  microcomputer  via  one  bit  of  a  parallel 
input/output  port.  Since  G(r,t)  is  available  to  the  micro¬ 
computer  as  a  logic  level  signal,  an  analog  to  digital 
converter  is  unnecessary  and  G(r,t)  is  merely  sensed  at 
intervals  controlled  by  the  external  clock. 

2.  The  Software 

Rules  #2  and  #3  of  Appendix  B  for  the  formulation  of 
I(r,t)  concern  the  application  of  the  hold  time  and  are 
repeated  below: 

2.  All  bursts  and  voids  must  have  a  duration  greater 
than  to  be  included  in  I(r,t). 

3.  For  the  conflicting  case  where  a  series  of  short 
bursts  are  separated  by  short  voids,  all  of  which  are 
less  than  t^,  a  single  uninterrupted  burst  will  result 
in  I(r,t)  provided  that  the  total  duration  is  then 
longer  than  t^. 

The  difficulty  in  applying  the  rules  is  that  they 
require  future  knowledge  of  the  fate  of  each  burst  and  void 
before  I(r,t)  can  be  formulated.  Rule  2  requires  the  ability 
to  "look  ahead"  msecs.  While  Rule  3  requires  as  much  as 
2t.  msecs,  of  future  knowledge  to  determine  the  instantaneous 
state  of  I (r , t) . 

The  "look  ahead"  capability  has  been  incorporated  into 
the  present  system  through  the  use  of  a  circular  data  buffer 
and  the  allowance  of  a  constant  time  lag  between  input  G(r,t) 
and  output  I(r,t  -  2 t ^ ) .  Figure  C-2  illustrates  the  concept 
of  a  circular  data  buffer.  The  three  data  pointers  are 
allowed  to  step  through  256  consecutive  bytes  of  memory 
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which  serve  as  a  data  buffer.  As  each  pointer  reaches  the 
end  of  the  buffer,  MEMORY(255),  it  is  automatically  reset 
to  MEMORY (0)  and  allowed  to  step  through  the  buffer  again 
in  a  circular  manner. 

At  each  clock  pulse,  new  G(r,t)  data  is  input  and 
stored  at  M(IPTR)  and  updated  I(r,t  -  2t^)  data  is  output 
from  M(OPTR) .  Following  this,  IPTR  and  OPTR  are 
incremented  so  that  they  are  ready  for  the  next  clock  pulse 
while  maintaining  the  same  relative  position.  The  output 
I(r,t  -  2x^)  always  lags  the  input  G(r,t)  by  a  constant 
2t.  msecs. 

“  -* 

Following  the  input/output  operations,  the  new  G(r,t) 

data  is  combined  with  previous  information  for  the  purpose 
of  measuring  burst  and  void  durations  and  comparing  them 

with  t.  .  If  sufficient  data  is  unavailable  to  determine  the 

h 

short  term  future  of  a  burst  or  void,  the  analysis  segment 
of  the  algorithm  pauses  and  waits  for  the  next  clock  pulse 
and  additional  G(r,t)  data.  However,  if  sufficient  infor¬ 
mation  is  available,  completely  formulated  I(r,t  -  nAt) 
data  is  inserted  into  the  buffer  starting  at  the  memory 
location  pointed  to  by  the  UPDATE  pointer.  Although  the 
UPDATE  pointer  is  always  bracketed  by  IPTR  and  OPTR,  its 
relative  position  fluctuates  depending  on  the  random  analy¬ 
sis  sequence. 

The  machine  language  routine  which  performs  the  hold 
time  application  is  given  at  the  end  of  this  section.  A 
unique  feature  of  the  RCA  1802  microprocessor  is  an  internal 
array  of  16  general  purpose  registers  which  may  be  used 
interchangeably  as  memory  pointers,  program  counters,  or  for 
temporary  storage.  Table  C-l  gives  the  register  assignments 
used  in  the  analysis  routines.  The  entry  point  for  the 
routine  is  at  label  "INITIAL"  or  (04  00)  and  it  is  in  this 
section  that  the  registers  are  initialized. 


\ 
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REGISTER 

0 

1 

2 

3 

4 

5 

6 
7 

8.1 
8.0 
9.1 
9.0 
A.  1 


Table  C-l  Register  Assignments 

ASSIGNMENT 


UPDATE  routine  program  counter 
Stack  pointer 

MAIN  routine  program  counter 
TAP  routine  program  counter 
IPTR  -  Input  buffer  pointer 
OPTR  -  Output  buffer  pointer 
UPTR  -  Update  buffer  pointer 
Constant  -  Buffer  page  number 
NB  -  Burst  counter 
Temporary  storage  for  I ( t) 

NV  -  Void  counter 
Constant  -  01 


A.  0 


Constant  -  00 
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The  analysis  program  is  made  up  of  a  main  routine, 

MAIN,  at  (0500)  and  two  subroutines,  TAP  (Turbulence 
Analysis  Program)  and  UPDATE  at  (0518)  and  (0570) , 
respectively.  MAIN  performs  the  functions  of  monitoring 
the  external  clock  and  input/output  control.  The  actual 
application  of  is  done  by  TAP  which  has  multiple  entry  and 
exit  points  depending  on  the  instantaneous  status  of  the 
analysis.  Once  the  bursts  and  voids  have  been  analyzed  by 
TAP,  UPDATE  performs  the  function  of  modifying  the  data 
buffer  in  accord  with  the  results.  For  additional  informa¬ 
tion  concerning  the  machine  language  programming,  the  reader 
is  referred  to  the  RCA  1802  Users  Manual. 

3.  Limitations 

The  precision  with  which  is  applied  to  G(r,t)  is 
limited  by  the  resolution  with  which  the  duration  of  the 
bursts  and  voids  may  be  determined.  Therefore,  it  is 
desirable  to  operate  the  external  sampling  clock  at  the 
highest  possible  frequency  compatible  with  the  data  handling 
capacity  of  the  microcomputer.  Each  possible  sequence  of 
instructions  that  the  microcomputer  could  be  required  to 
perform  during  a  single  sampling  interval  AT  was  examined 
and  the  time  required  for  their  execution  was  determined. 

The  maximum  sampling  frequency  could  then  be  chosen 
according  to  the  longest  execution  time  requirement.  It  is 
possible  to  increase  the  timing  resolution  by  increasing  the 
microcomputer  machine  cycle  speed,  thereby  improving  its  data 
handling  capacity  and  allowing  a  higher  sampling  frequency. 
The  microcomputer  utilized  in  the  present  experiment  was 
operated  at  a  machine  cycle  frequency  of  .25  Mhz. 

Table  C-2  lists  the  values  for  x.  and  the  sampling 

**  -V 

interval  AT  used  in  making  the  measurements  of  y(r)  and 
F^(r).  The  ratio  AT/x^  provides  a  measure  of  the  relative 
precision  with  which  the  hold  time  could  be  applied.  Due  to 
the  speed  limitations  of  the  microcomputer,  this  precision 
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Table  C-2  Signal  Analysis  Timing  Parameters 


x/D 

msec 

AT  msec 

AT/Th 

20 

2.15 

.54 

.25 

30 

4.0 

.67 

.17 

40 

6.19 

.77 

.12 

50 

8.68 

.87 

.10 

60 

11.44 

.95 

.08 

) 


becomes  progressively  more  critical  at  positions  closer  to 
the  jet  exit  and  accounts  for  the  discrepancies  evident  in 
Figure  C-3.  Accordingly,  the  hybrid  instrument  was  only 
used  for  measurements  at  x/D  >  20  . 

4 .  A  Microcomputer  Delay  Line 

For  correlations  involving  intermittency  functions 
with  relative  time  shifts,  a  second  RCA  based  microcomputer 
was  used  to  introduce  the  appropriate  delay. 

The  input  intermittency  signal  is  sampled  at  5  kHz 
rate  with  the  data  going  into  a  circular  data  buffer  at  a 
location  pointed  to  by  the  input  pointer.  At  the  same  time, 
data  from  the  buffer  pointed  to  by  the  output  pointer  is 
output.  As  the  input  and  output  pointer  are  stepped 
through  the  buffer,  the  output  signal  is  an  image  of  the 
input  but  delayed  by  an  amount  determined  from  the  separation 
between  the  pointers.  The  machine  language  program  which 
accomplishes  this  function  is  given  at  the  end  of  this 
section. 


t 
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Table  C-3 

TURBULENCE  ANALYSIS  ROUTINES 


ADDRESS 

CODE 

LABEL 

MNEM. 

OPERAND 

COMMENTS 

O'*  00 

F8  04 

INITIAL 

LDI 

04 

02 

B3 

PHI 

3 

03 

81 

PHI 

1 

04 

F8  06 

LDI 

PC  3 

06 

A3 

PC  3 

PLO 

3 

R3  =  Program  Ctr 

07 

03 

SEP 

3 

08 

F8  03 

LDI 

03 

OA 

B2 

PHI 

2 

OB 

f8  FF 

LDI 

FF 

OD 

A2 

PLO 

2 

OE 

E2 

SEX 

2 

R2  =  Stack  Ptr 

OF 

6C 

INP 

4 

Input  NHT 

10 

64 

OUT 

4 

Display  NHT 

11 

22 

DEC 

2 

M (RX)  -  NHT 

12 

F8  31 

LDI 

31 

14 

A1 

PLO 

1 

Rl  =  UPDATE  Pro¬ 
gram  Ctr. 

15 

F8  05 

LDI 

05 

17 

B4 

PHI 

4 

18 

F8  19 

LDI 

19 

1A 

A4 

PLO 

4 

R4  =  TAP  Pro- 

gram  Ctr. 


152 


Table  C-3  cont. 


ADDRESS 

CODE 

LABEL 

MNEH. 

OPERAND 

COMMENTS 

04  IB 

F8  06 

LDI 

06 

ID 

B5 

B5 

IE 

B6 

B6 

IF 

B7 

PHI 

7 

20 

B8 

PHI 

8 

21 

F8  30 

LDI 

30 

23 

A5 

PLO 

5 

R5  =  IPTR 

24 

A7 

PLO 

7 

R7  =  UPTR 

25 

F8  00 

LDI 

00 

27 

A6 

PLO 

6 

R6  -  OPTR 

28 

A8 

PLO 

8 

R8.0  =  NB 

29 

B9 

PHI 

9 

2A 

A9 

PLO 

9 

R9.0  =  NV 

2B 

AA 

PLO 

A 

RA.O  =  Con¬ 
stant  = 

2C 

F8  01 

LDI 

01 

2E 

BA 

PHI 

A 

RA .  1  =»  Con¬ 
stant  ® 

2F 

CO  0<*00 

LBR 

MAIN 

Jump  to  MAIN 

05  00 

3E  00 

MAIN 

BN3 

MAIN 

Wait  for  clock 
pulse 

02 

E6 

SEX 

6 

03 

61 

0UT1 

Output  M(OPTR) 

OPTR  +  1 
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Table  C-3  cont. 

ADDRESS 

CODE 

LABEL  MNEM. 

OPERAND 

COMMENTS 

05  0*» 

E5 

SEX 

5 

05 

69 

INN 

Input  data  to 
M(IPTR) 

06 

E2 

SEX 

2 

07 

Ok 

SEP 

k 

Jimp  to  TAP; 

D  =  G(t) 

08 

15 

INC 

5 

IPTR  +  1 

09 

98 

GHI 

8 

0A 

85 

PHI 

5 

OB 

B6 

PHI 

6 

OC 

36  OC 

WAIT  B3 

WAIT 

Wait  until  end 
of  clock  cycle, 
then  repeat 

OE 

30  00 

BR 

MAIN 

05  18 

03 

B1 

SEP 

3 

Return  to  MAIN 

19 

32  30 

NTRB1 

BZ 

V2 

If  G (t)  =  0,  go 
to  V2,  otherwise, 
NB  +  1 

IB 

18 

B2 

INC 

8 

1C 

88 

GLO 

8 

ID 

F7 

SM 

IE 

3B  18 

BL 

Bl 

If  NB  <  NHT,  re¬ 
peat 

20 

9A 

GHI 

A 

21 

89 

PHI 

9 

01  -►  1  for  UP- 

DATE 
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ADDRESS 

CODE 

Table  C-3  cont. 
LABEL  HNEM. 

OPERAND 

COMMENTS 

05  22 

88 

GLO 

8 

NB  D 

23 

D1 

SEP 

1 

Call  UPDATE 

24 

8a 

GLO 

A 

25 

A8 

PLO 

8 

0  ->  NB 

26 

D3 

B3  SEP 

3 

Return  to  MAIN 

27 

32  30 

NTRB2  BZ 

V  2 

If  G (t)  =  0,  go 
to  V 2,  otherwise, 
01  •>  D 

29 

9A 

GHI 

A 

2A 

D1 

SEP 

1 

Call  UPDATE,  then 
repeat 

26 

30  26 

BR 

B3 

2D 

D3 

VI 

SEP 

3 

Return  to  MAIN 

2E 

3A  ue 

NTRVl 

BNZ 

B4 

If  G(t)  *  0,  go 
to  B*»,  otherwise, 
NB  +  1  and  NV  +  1 

30 

18 

V2 

INC 

8 

31 

19 

INC 

9 

32 

89 

GLO 

9 

33 

F7 

SM 

3«* 

3B  2D 

BL 

VI 

If  NV  <  NHT,  re¬ 
peat 

36 

8a 

GLO 

A 

37 

B9 

PHI 

9 

00  -*  1  for  UPDATE 

38 

88 

GLO 

8 

NB  D 
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Table  C-3  cont. 


ADDRESS 

CODE 

LABEL 

MNEM. 

OPERAND 

COMMENTS 

05  39 

D1 

SEP 

1 

Call  UPDATE 

3A 

8A 

GLO 

A 

CD 

a8 

PLO 

8 

0  +  NB 

3C 

A9 

PLO 

9 

0  -*■  NV 

30 

D3 

V3 

SEP 

3 

Return  to  MAIN 

3E 

3A  IB 

NTRV2 

BNZ 

B2 

If  G(t),  to  to  82 

1*0 

9A 

GHI 

A 

NV  =  01  -*■  D 

1*1 

D1 

SEP 

1 

Call  UPDATE,  then 
repeat 

1*2 

30  3D 

BR 

V  3 

05  1*6 

18 

B4 

INC 

8 

NB  +  1 

47 

8a 

GLO 

A 

1*8 

A9 

PLO 

9 

0  -*•  NV 

49 

07 

LDN 

7 

If  M(UPTR)  -  1, 
go  to  B5 

4A 

3A  5A 

BNZ 

B5 

4C 

88 

GLO 

8 

1*0 

F7 

SM 

l*E 

3B  18 

BL 

B1 

If  NB  <  NV,  go  to 
81 ,  otherwise. 

50 

9A 

Gil’ 

A 

51 

B9 

PHI 

9 

01  -*•  1  for  UPDATE 

52 

88 

GLO 

8 

NB  ■*>  D 

156 


Table  C-3  cont. 


ADDRESS 

CODE 

LABEL 

MNEM. 

OPERAND 

COMMENTS 

05  53 

Dl 

SEP 

1 

Cali  UPDATE 

5A 

8A 

GLO 

A 

55 

A8 

PLO 

8 

0  -+■  NB 

56 

30  26 

BR 

B3 

then  go  to  B3 

5A 

9A 

B5  GHI 

A 

5B 

B9 

PHI 

9 

01  +  1  for  UPDATE 

5C 

88 

GLO 

8 

NB  -*•  D 

5D 

Dl 

SEP 

1 

Call  UPDATE 

5E 

8a 

GLO 

A 

5F 

A8 

PLO 

8 

0  NB 

60 

30  26 

BR 

B3 

then  go  to  B3 

05  70 

DA 

UEXIT 

SEP 

1» 

71 

22 

NTRU 

DEC 

2 

72 

52 

STR 

2 

73 

FE 

SHL 

7A 

FE 

SHL 

75 

FA 

ADD 

76 

FD  BA 

SOI 

10 

78 

A1 

PLO 

1 

Return  to  TAP 
RX  -  1,  Save  NHT 
NB  or  NV  •+  M(RX) 

Compute  address 
for  correct  number 
of  UPDATE  opera¬ 
tions 


Jump  to  address 
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Table 

C-3 

cont. 

ADDRESS 

CODE 

LABEL 

MNEM 

OPERAND 

05  79 

17 

U13 

INC 

7 

7A 

98 

GHI 

8 

7B 

B7 

PHI 

7 

7C 

99 

GHI 

9 

7D 

57 

STR 

7 

7E 

17 

UI2 

INC 

7 

etc. 

• 

05  B5 

17 

Ul 

• 

INC 

7 

B6 

98 

GHI 

8 

B7 

B7 

PHI 

7 

B8 

99 

GHI 

9 

B9 

57 

STR 

7 

BA 

12 

UO 

INC 

2 

BB 

30  70 

BR 

UEXIT 

COMMENTS 

UPTR  +  1 

I  -*■  M(UPTR) 


Sequence  repeats 
for  U 1 2  -  Ul 


UPTR  +  1 

I  *  m(uptr) 

Restore  M(RX)  * 
NHT 

Go  to  UEXIT 
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Table  C-4 
DELAY  LINE  ROUTINE 


ADDRESS 

CODE 

LABEL 

MNEH. 

OPERAND 

COMMENTS 

10  00 

F8  00 

START 

LDI 

00 

02 

B4 

PHI 

4 

03 

A4 

PLO 

4 

04 

E4 

SEX 

4 

R4  -  OPTR 

05 

7B 

SEQ 

Signal  for  input 

06 

3F  06 

A1 

BN4 

A1 

08 

6C 

1NP 

C 

09 

FA  OF 

ANI 

OF 

OB 

54 

STR 

4 

OC 

B5 

PHI 

5 

Input  and  display 

OD 

64 

OUT 

4 

High  delay  byte 

OE 

24 

DEC 

4 

OF 

37  OF 

A2 

B4 

A2 

11 

37  11 

A3 

BN4 

A3 

13 

6C 

INP 

C 

R5  **  IPTR 

14 

64 

OUT 

4 

Input  and  display 

16 

A5 

PLO 

5 

Low  delay  byte 

17 

7A 

REQ 

18 

37  18 

A4 

B4 

A4 

1A 

3E  1A 

LOOP 

BN3 

LOOP 

Wait  for  clock 
pulse 
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Table  C-4  cont. 


ADDRESS 

CODE 

LABEL 

MNEM. 

OPERAND 

10  1C 

E5 

SEX 

5 

10 

6f 

INP 

F 

IE 

64 

OUT 

4 

IF 

E4 

SEX 

4 

20 

67 

OUT 

7 

21 

94 

GHI 

4 

22 

FA  OF 

ANI 

OF 

24 

84 

PHI 

4 

25 

95 

GHI 

5 

26 

FA  OF 

ANI 

OF 

28 

B5 

PHI 

5 

29 

3F  1A 

BN4 

LOOP 

2B 

30  00 

BR 

START 

COMMENTS 


Input  and  display 
Current  data 


Output  delayed 
data 


Repeat  unless 
delay  change 
desired 


APPENDIX  D 


NUMERICAL  SENSITIVITY  STUDY 

In  the  formulation  of  the  numerically  simulated  vortex 
street  structure  given  in  Chapter  V,  it  was  necessary  to 
choose  values  for  four  similarity  proportionality  constants 
and  to  select  a  model  for  the  vortices.  The  effects  of 
variations  from  the  choices  which  were  utilized  in  Chapter  V 
are  given  in  this  appendix. 

1.  Proportionality  Constants 

Of  the  four  proportionality  constants  involved  in  the 
numerical  simulation,  the  pattern  wavelength  coefficient, 

CL#  and  the  vortex  strength  coefficient,  Cr ,  could  be  set 
directly  from  experimental  observations.  The  choices  for 
the  offset  distance  coefficient,  CH,  and  vortex  core  radius 
coefficient,  CR,  were  based  on  less  objective  criteria. 

For  these  two  parameters,  the  selection  was  based  upon 
finding  the  optimum  combination  which  would  yield  the  best 
agreement  between  numerically  calculated  and  experimentally 
determined  distributions  of  U,  V  and  UV.  The  values 
chosen  were  CH  =  ,05  and  CR  =  .13. 

To  examine  the  sensitivity  of  the  calculations  to 
these  parameters,  they  were  repeated  with  +  25%  variations 
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1  1 


of  the  coefficients  from  their  "optimum"  values.  Only  one 
coefficient  was  varied  for  each  case. 

Figures  D-l  through  D-4  illustrate  the  sensitivity  of 
the  calculations  to  CR  and  Figures  D-5  through  D-8  demon¬ 
strate  the  effects  of  variation  of  C„.  From  the  figures,  it 
is  noted  that  CR  has  a  comparatively  stronger  influence  on 
the  mean  velocity  profiles  while  the  Reynolds  stress  and 
fluctuating  intensities  are  more  sensitive  to  CH. 

2 .  The  Vortex  Model 

For  the  calculations  described  in  Chapter  V,  a  Rankine 
vortex  was  utilized  to  represent  the  individual  vortices 
making  up  the  vortex  street.  In  the  discussion  of  the 
results,  it  was  argued  that  the  capability  of  the  model  to 
calculate  reasonable  distributions  of  mean  lateral  velocity 
and  Reynolds  stress  is  a  direct  consequence  of  the  model 
being  able  to  calculate  the  mean  longitudinal  velocity 
accurately. 

Figure  D-9  illustrates  the  tangential  velocity 
distribution  for  an  alternate  choice  of  the  vortex  model. 

In  this  representation,  the  induced  velocity  is  exactly  a 
potential  vortex  for  r  >  rc  and  solid  body  rotation  for 
r  <  rc.  The  discontinuity  of  the  velocity  derivative  at 
r  =  rc  for  the  alternate  choice  should  result  in  significant 
changes  in  the  calculated  results  if  it  hampers  the 
capability  of  the  model  to  predict  the  mean  longitudinal 
velocity. 


Figure  D-6  Sensitivity  of  V(X,Y/B)/U  (X)  to 


Figure  D-7  Sensitivity  of  |UV  (X, Y/B) /U 
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The  results  of  calculations  utilizing  the  alternate 
vortex  are  given  in  Figures  D-10  through  D-13.  For  these 
calculations ,  and  CR  were  held  to  the  same  values  as  for 
those  with  the  Rankine  vortex.  Cp  was  set  to  2.401  rather 
than  2.597  in  order  to  match  the  experimentally  determined 
mean  centerline  velocity.  A  value  of  CR  =  .125  instead  of 
CR  =  .130  was  needed  to  provide  the  expected  spreading  rate 
(i.e. ,  U/tfm  =  .5  at  Y/B  =  1). 

From  the  results,  it  is  noted  that  the  subtle  changes 
in  the  slope  of  the  mean  longitudinal  velocity  profile  are 
associated  with  changes  in  the  lateral  velocity  and 
Reynolds  stress.  The  overall  magnitudes  are,  however, 
approximately  the  same  as  for  the  calculations  using  Rankine 
vortices.  The  "saddle-back"  distribution  of  the  Reynolds 
stress  might  have  been  expected  from  the  presence  of  the 
velocity  derivative  term  in  Eq.  5.57. 


s,n*s” 


Figure  D- 12  Sensitivity  of  |UV(X,Y/B 

Vortex  Model  Choice 
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